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INTRODUCTION 


Considerable work has been performed over the past 15 years on the 
development of orthotropic photoelastic material and the development of 
stress-optic laws which characterize the relations between the optical 
and mechanical responses of the material. Some plane stress elasticity 
problems have been solved with the material but these problems have been 
more or less of a demonstration nature, designed to show the extent to 
which orthotropic photoelasticity works. In many cases the problems 
solved were such that the load was applied parallel to or perpendicular 
to the fiber direction, the principal stress system coincided with the 
principal material system, the principal stress directions were known, 
the stresses were computed along paths parallel or perpendicular to the 
fiber direction, or the vast majority of the boundaries were traction 
free so that 2 out of the 3 stresses were known with certainty over much 
of the boundary. Oftentimes the problems were chosen because there ex- 
isted known solutions, or good first approximations, to the state of 
stress in the body. The stresses calculated using orthotropic photo- 
elasticity were compared to the known solution. These types of exerci- 
ses were necessary in the development stages of this new technology. 
However, to have the same utility as classical isotropic photoelastici- 
ty, orthotropic photoelasticity must be applicable to a wide variety of 
situations. Ultimately, orthotropic photoelasticity must be able to pre- 
dict the state of stress when the applied load is at an arbitrary angle 
relative to the fibers, when the principal stress directions are vari- 
able from point to point and are unknown, and when the stresses along 
arbitrary paths are required. In addition, if residual birefringence 
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exists in the material due to the manufacturing process, its effect on 
interpreting photoelastic data must be taken into account. 

This study was initiated to assess the ability of currently accept- 
ed orthotropic stress-optic laws to accurately predict the state of 
plane stress in currently available orthotropic photoelastic material. 
The ability to be able to examine a variety of loading directions, rela- 
tive to the axes of orthotropicity, was of prime concern. In addition, 
it was important to be able to compute the stresses at more or less ar- 
bitrary locations in the interior of the body. In classical isotropic 
photoelasticity, if one of the principal stresses is known, the other 
principal stress can be computed directly from the isochromatic fringe 
pattern. This is always the case when examining the stresses at a trac- 
tion-free boundary. For this case the shear and normal stresses are 
zero. With the shear stress zero, the principal stress directions are 
tangent and normal to the traction-free edge. With the normal stress 
zero, the only unknown is the principal stress tangent to the free edge. 
Using the stress-optic coefficient of the material and the isochromatic 
fringe count, this non-zero principal stress can be computed. In the 
interior of a body, both the isochromatic and isoclinic fringe patterns 
are used to compute stresses. Unless a third point-by-point experiment- 
al measurement is used in conjunction with the two photoelastic fringe 
measurements, an auxiliary numerical scheme is generally used to compute 
the complete state of stress in the interior of the body. This auxili- 
ary numerical scheme is usually in the form of the numerical integration 
of the plane-stress equilibrium equations or the compatibility equations 
along a path from a boundary, or other region of known conditions, to 
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the interior point of interest. For simplicity, this path is usually a 
straight line. 

This report begins by developing the stress-optic law for orthotro- 
pic materials. The laws are developed on a physical basis. Some of the 
development draws from the past work in orthotropic photoelasticity and 
can be considered as a summary and/or a reformating of results. How- 
ever, early in this study it became obvious that the techniques used to 
fabricate the photoelastic material could produce a residual birefrin- 
gence. With no applied loads, there could be an isochromatic fringe 
count of, say, 0.1. There would, of course, be an isoclinic associated 
with this value. The development which follows includes the effect of 
residual birefringence in the stress-optic law. Conceptually, account- 
ing for residual bi refringence effects is simple. Accounting for resid- 
ual effects operationally is also quite simple. Since bi refringence is 
associated with tensor quantities, keeping track of residual effects is 
bookkeeping. 

To have the ability to generate a wide variety of stress states ex- 
perimentally and at the same time to have a good approximation to the 
state of stress, an unidirectional disk in diametral compression was 
chosen as the study specimen. With such a specimen the fibers could be 
oriented at arbitrary directions relative to the applied diametral load 
and a wide variety of stress states could be produced. In addition, the 
stress state along arbitrary diameters or chords could be computed using 
the photoelastic data and an auxiliary numerical scheme. It was decided 
to refrain from using a third point-by-point experimental measure to 
compute the stresses. The most accurate point-by-point measure is to 
use interferometry to measure the thickness change of the material. 
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This requires optical equipment beyond the standard polari scope set-up. 
It was a goal to keep this study within the context of common polari- 
scopes. Oblique incidence techniques were not considered because of the 
poor fringe resolution associated with orthotropic photoelastic materi- 
als. 

The disk geometry was chosen for several other reasons. First, it 
was an easy specimen geometry to make. Second, and more importantly, it 
did not use much of the valuable photoelastic material that was avail- 
able. One specimen could be used for all load orientations. Finally, 
the disk geometry is an ideal specimen for calibration of orthotropic 
photoelastic material. In theory, according to the state-of-the-art 
stress-optic laws, three tests must be made on an orthotropic photoelas- 
tic material to completly characterize its stress-optic relations. In 
the past, three tensile tests have been conducted. One tensile test is 
conducted on a 0° specimen, one test on a 90° specimen, and one test on 
a 45° specimen. The experimenter can judiciously choose the length and 
width of the tensile specimens but the specimens do tend to consume the 
valuable photoelastic material. Since the material has some variability 
in its mechanical and optical properties, it is a good idea to calibrate 
the actual panel from which experimental models are to be made. Al- 
though three is a minimum, conceivably more tensile specimens, say at 
30°, 60°, and 75°, could be tested and the calibration could be done 
with some form of least-squares approach to the data reduction. This 
would be done, of course, at the expense of valuable photoelastic mate- 
rial. This has not been done, rather, the minimum three tests have gen- 
erally been the standard. With a single disk, theoretically, an infi- 
nite number of stress states is possible. Practically speaking, with 
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various disk orientations relative to the load direction, and with pho- 
toelastic data taken at several points on the disk, 15-20 sets of stress 
and corresponding photoelastic fringe information could be readily ob- 
tained and a least-squares calibration could be conducted to determine 
the stress-optic characterization of the material. 

Experimentally there were many compelling reasons for choosing the 
disk geometry. In addition, the theoretical elasticity solution to the 
disk in diametral compression could be derived from the complex variable 
approach. The complex variable solution is presented in an appendix and 
will not be discussed further except as to quote stress predictions from 
the solution. 

This report summarizes the work performed using the disk as a study 
specimen. The work followed the basic line of reasoning: If a particu- 

lar stress-optic law for the orthotropic photoelastic material is true, 
and if the calibration of the material is done using three tensile 
tests, then two types of experiments can be done to test the validity of 
the stress-optic law. Using the calibration information, the stress- 
optic law, and the theoretical elasticity solution, the isochromatic and 
isoclinic fringe patterns in the disk can be predicted. By examining 
different load orientations, relative to the fiber directions, a wide 
variety of spatial variations of fringe patterns can be generated. This 
wide variety of predicted fringe patterns can then be compared to the 
experimentally observed fringe patterns. This is a legitimate test of 
the assumed stress-optic behavior and the calibration. These so-called 
"forward tests" were conducted and are reported. Real interest, of 
course, is in starting with a specific fringe pattern and working back- 
wards to obtain the stress state. This is the essence of photoelastici- 
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ty. A numerical scheme was previously developed for studying the stres- 
ses in the inner laps of double-lapped, double-pin connector models con- 
structed from acrylic and isotropic photoelastic material. This scheme 
was based on the stress-optic law for classic isotropic materials and a 
finite-difference form of the plane-stress equilibrium equations. The 
equilibrium equations were chosen as the auxiliary conditions to obtain 
the three components of stress at interior points in the model. This 
scheme was modified to accomodate an orthotropic stress-optic law. The 
stress-optic calibration information from auxiliary calibration tests, 
the isochromatic and isoclinic fringe count at the various mesh points 
of the finite-difference grid, and the finite-difference equilibrium 
equations were used to experimentally determine the stresses in the 
disk. These stresses were compared with the stresses computed from the 
complex variable elasticity solution. These stress prediction tests 
could be considered as "backwards tests." To put the results into con- 
text, numerical and experimental results for disks made from isotropic 
photoelastic material are compared. Also numerical studies were conduc- 
ted to determine the effect of isocline measurement error on the experi- 
mental results. 

Throughout the report specific references to previous investiga- 
tions are made. These are noted in a list of references. However, for 
completeness, a bibliography is included at the end of this report to 
mention those investigations not specifically cited. 
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ORTHOTROPIC PHOTOELASTICITY 


Bi refringence 

In many transparent solids, the index of refraction* can be consid- 
ered as a quantity which is independent of the direction of propagation 
and the plane of polarization of the transmitted light. However, in 
some transparent solids the speed of propagation of light depends both 
on the direction of propagation and the plane of polarization of the 
light. Within such a solid there will be a direction of propagation and 
a plane of polarization for which the light travels the fastest. There 
is also a direction of propagation and plane of polarization for which 
light travels the slowest. Because of the orthogonality and symmetry of 
the molecular structure and crystals making up these solids, many of 
these directions and planes are perpendicular to each other. This idea 
of multiple propagation speeds and directional -dependence of the speed 
of propagation is termed bi refringence . In the discussions here, the 
birefringence phenomenon will be considered something that can occur in 
an unstressed material and also something that can be induced in a ma- 
terial by application of loads. Because of the orthogonality of the 
directions associated with the different propagation speeds, the indices 
of refraction can be associated with quantities which can be regarded as 
tensors. Tensor quantities have principal values and the directions as- 
sociated with these values are quite often perpendicular to each other. 
The ideas of extremums in the propagation speeds, tensors, planes of po- 
larization, etc., can all be used to put the ideas of classical 


*index of refraction of a medium, n, is, defined as the speed of light in 
a vacuum, c Q , divided by the speed of light in the medium, c. 
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isotropic photoelasticity on a formal and rigorous foundation. This has 
been done in various texts and papers. However, if one is trying to de- 
termine, for example, the stress concentration factor at the base of a 
fillet in a model representing a metal component, the formalism is not 
required. All that is required is to know that the isochromatic fringe 
number times some a priori determined numerical factor gives the value 
of the stress at the base of the fillet. However, if one is dealing 
with a new photoelastic phenomenon such as orthotropic photoelasticity, 
the formalism helps guide the interpretation of the photoelastic data. 
The formalism also provides insight into the potential of the new phe- 
nomenon. In this report, some of the formalism will be used to estab- 
lish ideas regarding orthotropic photoelasticity. The dielectric tensor 
and its relation to the indices of refraction will be introduced. Then 
the notion of the bi ref ringence tensor will be discussed. The introduc- 
tion of the birefringence tensor is a convenient abstraction and is not 
at all necessary in order to understand the implications of orthotropic 
photoelasticity. However, many past investigators have used the idea of 
a birefringence tensor and so it will be used here to illustrate how it 
fits into the nomenclature and mechanics of orthotropic photoelasticity. 
These ideas will be used to derive relations between the applied stress 
and the observed optical response. These relations will be such that 
they can be used in a numerical scheme to compute the stresses in a 
loaded photoelastic model. To begin, however, the approach used by many 
texts to start the discussion of isotropic photoelasticity will be re- 
produced. The expression for the relative retardation will be derived. 
This approach will introduce the relations between the wavelength of 
light, the model thickness, the indices of refraction, and the slowing 
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of the light passing through photoelastic material. For those not 
familiar with even isotropic photoelasticity, this approach will be 
quite palatable because it is based on the ideas of travel times of a 
light wave, optical path lengths, etc. Also, since many texts use this 
idea, it will be a point of reference if one is consulting these texts. 
For those familiar with isotropic photoelasticity, the derivation will 
be a review and it will provide some commonality between isotropic and 
orthotropic photoelasticity. The relative retardation is introduced be- 
cause it is one of the two basic quantities measured directly by a 
standard polar i scope. The other quantity, which will be introduced 
after the relative retardation is discussed, is the isocline. 

Relative Retardation 

Consider a transparent solid in the form of a thin plate, of thick- 
ness h, oriented perpendicular to the direction of propagation of a 
light vector. As shown in fig. 1, a light vector is impinging on or in- 
cident to the plate. If the light vector is polarized, it can be 
thought of as oscillating harmonically in a single plane as it propa- 
gates toward the plate, as shown in fig. 2. If the plate is made of a 
material with indices of refraction which are dependent on the plane of 
polarization of the incident light, then it is possible to identify a 
plane of polarization in which the light travels through the plate with 
the fastest propagation velocity. It is also possible to identify a 
plane of polarization in which the light passes through the plate with 
the slowest propagation velocity. These directions, as shown in fig. 3, 
are referred to as fast and slow directions, respectively. The descrip- 
tions fast and slow can be misleading. In a transparent solid light 
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Fig. 2 


Orientation of polarized light. 
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actually slows, when compared to its speed in air, independent of the 
plane of polarization. However, in the so-called fast direction the 
light slows the least. 

The polarized light that is impingent on the plate in some direc- 
tion other than the fast or slow direction, as shown in figure 4, can be 
decomposed into vector components in the fast direction and slow direc- 
tion, respecti vely. Referring to in fig. 5, due to differences in the 
speed of propagation in the two directions, these two components are 
out-of-phase with each other when they emerge from the plate. Since the 
components are out of phase they can be made to interfere constructively 
and destructively. The specific conditions required to obtain inter- 
ference is a function of the optical apparatus used for viewing the phe- 
nomenon. The various optical devices, or polari scope arrangements, and 
the specific relations for interference of these out-of-phase components 
are discussed in [1], It is this phase difference which is the basis 
for using the birefringent phenomenon to an advantage. An expression 
for this phase difference is developed in the next section. 

Phase Difference 

If the index of refraction in the fast direction is denoted as n 2 , 

then 



where c 0 is the speed of light in a vacuum and is the speed of light 
in the plate if the polarized light vector is aligned with the fast, 2, 
direction. The index of refraction in the slow, 1, direction is simi- 
larly defined, i.e. 
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Fig. 4 Decomposition of polarized light. 
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(2) 




The hat, A , over the indices of refraction signifies that these are the 
extremum values for all orientations of the plane of polarization. This 
hat notation will soon be used to denote the principal values of a ten- 
sor associated with the transmission of light through the plate. The 
principal values of this tensor and the extremums in the values of the 
indices of refraction are related. 

For light propagating with the fast speed, to pass through thick- 
ness h, the elapsed time, t 2 , is 


t 


2 



(3) 


Similarly, for light propagating with the slow speed, the elapsed time, 
t^, to pass through thickness h is 


t 



(4) 


If a vacuum occupied that region of space instead of a birefringent 
medium, then the elapsed time to move through the distance h would be 


t 


o 



(5) 


Figure 6 illustrates this thickness geometry. The time difference be- 
tween when the fast light emerges from the plate and when the slow light 
emerges from the plate. At, is given by 
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Fig. 6 Light passing through plate. 
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( 6 ) 


At = - tg 


h h 



Using eqs, 1-4 In eq. 6 leads to 

At = "r — (n, - n«) • (7) 

o 

In terms of distance (e.g. meters, inches, etc.), the two light vectors 
are out of phase by an amount 

c Q At = h (nj - n 2 ) . (8) 


Relative to the wavelength of light , X, entering the plate, the phase 
difference, In cycles, between the two light vectors is 


c 0 At 

X 




(9) 


To put this phase difference in terms of radians, multiply eq. (9) 
by 0 it, l.e. 


^ c o at ath 
X — (n l 



( 10 ) 


Thus the phase difference. A, is expressed in terms of the wavelength of 
Incident light, the plate thickness, and the difference in the principal 
indices of refraction. The quantity A is the relative retardation be- 
tween the two components of reemerging polarized light. 

A 

In this discussion the quantity n^ is arbitrarily associated with 
the direction that transmitts the light the fastest and n 1 is associated 
with the direction that transmitts the light the slowest. Thus 

A A 

(n^ - n 2 ) Is a positive quantity. 
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Causes of Birefringence in Solids 

The interaction of the electromagnetic field of the light with the 
electromagnetic fields associated with the bonding of the atoms in a 
transparent solid changes the speed of light from the vacuum condition. 
In addition, anything that changes the spatial distribution and/or 
strength of the bonding electric fields in a solid will affect the pro- 
pagation speed of light within the solid. A way to change the bonding 
electric fields in a solid is to physically move the molecules of the 
solid relative to each other. Such a condition could be caused by ap- 
plying a mechanical load (stress) to the solid. If the solid is elas- 
tically isotropic, and if there is a stress applied in a particular 
direction, then the strength and distribution of the bonding fields 
would be different in different directions. In this case the electric 
field of the incident light would be acted upon by one field strength if 
it was polarized in one direction and it would be acted upon by another 
field strength if it was polarized in another direction. This leads to 
di recti onally-dependent propagation speeds. This phenomenon would also 
occur if the material was elastically orthotropic. In fact, even with 
the solid in its natural (unstrained or unstressed) state, an orthotro- 
pic solid has different bonding field strengths in different directions. 
Thus it is possible to conceive of multiple indices of refraction in an 
orthotropic media with the molecules In their natural (unstressed) posi- 
tions. In the context of these discussions, the multiple Indices of re- 
fraction with no load applied can be considered a residual difference in 
the indices of refraction or residual birefringence. If a mechanically 
induced alteration of the bonds occurs, and if superposition of bire- 
fringent effects is assumed valid, then the total differences in indices 
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of refraction would be the sum of the effects of the residual bire- 


frlngence and the effects of the mechanically induced birefringence. 

Details of Birefringence in Solids 

The fundamental quantities for the study of birefringence In solids 
are the dielectric constants of the material. The plural connotation in 
the word 'constants' Is important because oftentimes the dielectric 
properties of a material are thought of as being described by a single 
number. In general, the dielectric properties of a material are direc- 
tional ly-dependent and they can be represented as a tensor called the 
dielectric tensor . The dielectric tensor relates the electric field 
strength vector to the electric displacement vector. The index of re- 
fraction of a solid is related to the dielectric tensor. Due to the di- 
rectionally-dependent nature of a tensor, the speed of propagation of 
light in a transparent solid depends strongly on the directional -depen- 
dence of the dielectric properties of the solid. For the case of a thin 
plate, referred to here by analogy to 'plane stress' as 'plane birefrin- 
gence,' the dielectric tensor can be written as 



With the above notation, the dielectric properties are being referred to 
an x-y coordinate system. The x-y coordinate system is a global coordi- 
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nate system with the x and y axes in the plane of the plate. At the mo- 
ment the orientation of the x and y axes within the plane are not impor- 
tant. Later, stresses will be computed in the x-y system and its orien- 
tation will be specified. Figure 7 depicts the coordinate system. The 
light can be considered to be propagating into the page and perpendicu- 
lar to it. 

In an electrically isotropic material the dielectric tensor takes 
the form 

e o 

o e 

A vacuum, or free space, is isotropic and so the dielectric tensor is 
given by 

e 
o 

the subscript 'o' denoting free vacuum. The fast and slow indices of 
refraction of a solid are related to the principal values of the 
dielectric tensor. The principal values of the dielectric tensor are 
given by 




e 2 = 


e x + e y 
— 2 


* 7 - 


e - e 
x 


*)2 + e 2 
J xy 


(14) 
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The fast and slow Indices of refraction are defined [2] to be 



(15) 

(16) 


For the transparent solids considered here and for the levels of 
birefringence encountered in photoelasticity, certain simplifying ap- 
proximations can be made. If e ave Is defined as 


then eq. 14 becomes 
where 


e x + e y 

e ave 2 ^ * 


A A 


e l» e 2 = e ave i Ae * 


(17) 

(IB) 


Ae = 




(19) 


It is shown in Appendix A that Ae is an order of magnitude smaller than 
e ave . Then, using the binomial approximation, eqs. 15 and 16 become 


A 


n 


1 * 


"2 



[1 ±le^"] * 


ave 


( 20 ) 

( 21 ) 


The relative retardation can be related to the dielectric constants 
by substituting eqs. 20 and 21 into eq. 10. This relation is important 
because the application of stress directly affects the dielectric con- 
stants of a transparent material. The above substitution results in 
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A 


( 22 ) 


- 2tt h Ae 
' Xn ave e o 


* 


where n ave 
defined as 


is the average index of refraction of the medium and is 


n 


ave 



(23) 


It should be mentioned at this point that it has been assumed that 
the magnetic permeability, y, has been assumed to be an isotropic prop- 
erty of the material. Furthermore, the permeability is assumed to be 
equal to the permeability in a vacuum, y Q . Generally the velocity of 
propagation of electromagnetic waves depends on the permeability as well 
as the dielectric constant. For an isotropic material, the velocity of 
propagation is proportional to l//ey . In a vacuum the velocity becomes 
l//e Q u 0 . The index of refraction of a medium is thus given by thus 



If y is assumed to be equal to y Q , then the index of refraction is 
related only by the ratios of the dielectric constants. 

For an orthotropic material, particularly fiber reinforced compos- 
ite materials, there is a coordinate associated with* the stiff (fiber) 
direction and a coordinate, perpendicular to the first, associated with 
the soft (matrix) direction. Much of the analysis of fiber-reinforced 
composites uses these directions as the reference for defining stress, 
strains, elastic constants, failure stresses, etc. In orthotropic pho- 
toelasticity it is also convenient to use this coordinate system as a 
basis. Herein, these directions will be referred to as 1M and 2M. 
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These directions will refer to the fiber and matrix direction, respec- 
tively, with the M denoting that it Is a material coordinate system* 

When computing stresses in the x-y system, the angle the +1M axis makes 
with to the +x axis will be denoted <J>* This coordinate system orienta- 
tion is shown in fig. 8. In the material coordinate system, the dielec- 
tric tensor can be identified as 



(25) 


Birefringence in Fiber-Reinforced Solids 

Because of differences in the thermoelastic properties of the fiber 
and matrix constituents in fiber-reinforced composites, the fabrication 
of fiber-reinforced composites can lead to residual stresses in the fi- 
bers and the matrix. Since the dielectric properties of solids can de- 
pend on the stress state in the solid, these residual stresses can lead 
to residual birefringence. As mentioned in the Introduction, such bire- 
fringence is considered here. In the material coordinate system, resid- 
ual dielectric effects generally have three components. 



(26) 


The superscript R identifies properties associated with residual ef- 
fects. There is no a priori reason to believe, in general, that the 


25 




principal directions of the residual dielectric tensor are aligned with 
the material directions. There Is, however, a principal direction asso- 
ciated with the residual tensor, i.e. directions associated with the 
fastest and slowest propagation directions for the solid when there Is 
no mechanical load applied. Assume that the principal residual direc- 
tions are rotated from the material directions by an amount e R . The or- 
ientation of the 1R-2R system relative to the 1M-2M system is shown in 
fig. 9. The residual dielectric tensor principal coordinate system Is 
identified by 1R and 2R. The principal values of the residual dielec- 

^ a. 

trie tensor are given by e 1R and e 2R . The representation of the residu- 
al dielectric tensor in the principal residual coordinate system is 

A I 

e lR I 

ej . (27) 

o I 


The components of the residual dielectric tensor in the material 
system, eq. 26, can be related to the principal values of the residual 
dielectric tensor, eq. 27, by the usual two-dimensional tensor trans- 
formations, i.e. 



e R 

e 12M 


r e lR + e 2R 
*■ 2 


r e lR + e 2R 

L 2 



.) ♦ t f A! 

0 - — ) 

— ) sin (-»„) 


COS (-29 r ) 
cos (-2e R ) 


(28)a 

(28)b 

(28)c 
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The minus sign is left with the angle as a reminder that in these equa- 
tions one is transforming back to the material system. If the compo- 
nents of the residual dielectric tensor are known in the material sys- 
tem, the principal values and principal directions of the residual diel- 
ectric tensor can be determined by the familiar formulas 


1M 


1R * 2R 


*2M 


A 


1M 


- e 


2M' 


( e 12M) 2 


(28)d 


2 

tan (20^) s ^ jj— • ( 28) e 

e lM e 2M 

If loads are applied to the solid and there are residual birefrin- 
gent effects present, then there will exist a state of _total_ birefrin- 
gence. This situation is associated with the total dielectric tensor. 
The total dielectric tensor will have principal values and principal 
directions. The coordinate system of the principal total dielectric 
tensor will be denoted as IT and 2T, the T signifying association with 
total effects. The principal values of the total dielectric tensor can 
be written as . 

I A 

yn 

/ A 

I e 2T 
\ 0 

In the material coordinate system, the components of the total dielec- 
tric tensor are given by 
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The angle between the material coordinates and the coordinate system of 
the principal total dielectric tensor is denoted as Qj and is illustra- 
ted in fig. 10. The components of the total dielectric tensor in the 
material coordinate system, eq. 30, are related to the principal values 
of the total dielectric tensor, eq. 29, by 


A A A 


e lM 

_ r IT e 2T>i w e lT " e 2T^ _ 

= { 2 J + 1 2 J cos 

A A A A 

(-2e x ) 

(31)a 

e T 

e 2M 

_ f e lT + e 2T-, r e lT ‘ e 2Ti 

l 2 j ( 2 J cos 

(-20-p) • 

( 31 ) b 

e T 

e 12M 

= - (-- - 2 - 62T ) sin (-20 t ) 


( 31 ) c 


Conversely, knowing the components of the total dielectric tensor in the 
material system* the principal values and principal directions of total 
effects can be computed by 


e lT’ e 2T 





( 31 ) d 


tan ( 20y) 




( 31)e 


Finally, there are principal values, principal directions, and ma- 
terial system components of the dielectric tensor associated with the 
application of mechanical loads. In the material system, the components 

of the dielectric tensor due to the applied loads are given by 
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Fig. 10 Orientation of 1T-2T system relative to the 1M-2M system. 
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The mechanically induced dielectric tensor has its own principal coordi- 
nate system 1 o - 2 a. The principal values of the dielectric tensor due 
to applied loads can be referred to this system and can be written as 



Assume the lo - 2 a coordinate system is oriented at an angle 0 relative 

a 

to the material system. Figure 11 shows the angle and four of the 
five coordinate systems discussed. (The x-y coordinate system is not de- 
picted.) Equations similar to eqs. 28a - 28e and 31a-31e can be written 
to relate e^, e^, e^, and e . These equations are 


e lM = ( 6lg 2 + ( 6 -— 2 — ) cos (-2e 0 ) { 34) a 

4 m = cos (-2e o ) ( 34 ) b 

A A 

e 12M = " 2 Sln { ' 29 a ) (34)C 





e, u + e. 


e° - e a 

Cm C/ 


„ * f 1M 2Mi j. \ If 1M 2 M ^9 f a i?a\a 

e lo* e 2a = C 2 J 1 VI 2 J + ( e 12MJ (34)d 

(34)e 


2e l?M 

tan (29 ) = 

O _0 _0 


e l« " e 2M 


Hypothesis Regarding Dielectric Effects Due to Applied Loads 

It will be hypothesised here that the dielectric effects due to 
applied loads can be superimposed onto the residual dielectric ef- 
fects. Due to the tensorial nature of dielectric effects, the total 
dielectric tensor can be represented, in the material system, as 
fol 1 ows : 


(35) 


/ T 
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[ R ) 

l e iM i 


1M 1 


e lM j 

) T ( 


) a \ 

| | a 

) R ( 

j e 2M l 


) e 2M 

' + < 

e 2M / 

It \ 

| j 
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„R \ 

\ e 12Mj 


, 12Mj 
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It has been hypothesized by previous investigators [3,4] that the 
dielectric tensor due to the applied loads is related to the applied 
stress by the following orthotropic relation: 


'1M 

} 2M 

12M, 


<11 <12 0 
<21 <22 0 
0 0 <66 


1 
°2 
° 12 , 


(36) 


The q-jj are constants for a given material and can be thought of as the 
material's stress-dielectric coefficients. Equation 35, representing 
the total dielectric tensor in the material system, can be written as 
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“11 

^12 

0 

j e 2M f = 

<21 

q 22 

0 

e T ) 

0 

0 

q 66_ 



Several approaches can be taken to check the validity of eq. 37 or 
to check the consequences of it, assuming it is valid. It is on the 
basis of eq, 36 that much of the work in orthotropic photoelasticity has 
been done to date. The hypothesis represented by eq. 36 has been termed 
ad hoc [5] because there is no a priori reason to believe that orthotro- 
pic elastic behavior implies orthotropic optical behavior. The crux of 
the discussion about the validity of eq. 36 is whether or not it corre- 
lates with experimental observation. Equation 35 has not been verified 
to any great extent. Previous investigators [6, 7, 8, 9] have mentioned 
residual birefringent effects. However eq. 35, and the just-derived 
relations involving the residual dielectric tensor, treat residual 
effects with more formalism. 


Consequences of the Hypothesized Relation 

Since, as stated previously, standard polariscopes measure relative 
retardation. A, it is useful to relate the retardation A to the applied 
stresses via the just stated hypotheses. This is done by using eq. 37 
in eq. 22, the expression for A. In the present context the dielectric 
properties in eq. 22 refer to the total dielectric properties, i.e. 
residual effects plus mechanically induced effects. Thus eq. 22 becomes 


A 


2irh 

T" 


Xn 


ave 



(38) 


35 



with 


and 


Ae T = 





(39)a 
( 39) b 
(39)c 


Since el w ^ is related to the trace of the total dielectric tensor, el.„ 
ave ave 

is independent of coordinate system. To use eq. 37 in eq, 39, it is 
useful to write the residual dielectric effects in eq, 37 in terms of 
principal residual values. This can be done using eqs. 28a - 28c. 

Doing this, eq. 37 becomes 



<11 

<21 

0 



A 



A A 

e lR‘ e 2R 


2—^) cos(- 29 r ) 

A 

cos(- 26 r ) 


e lR~ e 2R 


-) sin(-26 R ) 


(40) 


Substituting eq. 40 into eq. 39 and then substituting those results into 
eq. 38 leads to an expression for A, namely 

A 'TXT V t(q ll ' q 21 )o l - (q 22 ' q 12 )o 2 + <*1R - *2R> cos(29 rH 2 
Xn ave e o V 


+ ^ 2q 66 T 12 + 


1R 


' e 2R^ 


sin(2e„)} 2 


(41) 
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With a darkfield polari scope arrangement, for example, fringes occur 
when 

-ip - 0, it, 2ir , 3ir ,. . • ,Nit (42) 

N an integer. 

Combining eqs. 41 and 42 gives the condition for Integer fringe orders 
in a darkfield polariscope due to applied loads and residual effects. 
Specifically, integer order fringes occur when 

-, n T ' q 21 )o l ' * q 22 ■ q 12^2 + (e lR ’ e 2R^ cos ( 29 r)} 2 

^ An ave e o v 

+ t 2q 66 T 12 + ( ®1R - *2R } sin(20 R» 2 " N T 

(43) 

Ny being the integer fringe order. The subscript T has been added to N 
to signify it is the fringe order associated with total birefringence, 
residual and mechanically induced. Practically speaking, Ny does not 
have to be an integer. Standard polari scopes are capable of displaying 
any fringe order, fractional or integer. In any case, the residual ef- 
fects and the applied load effects combine In the manner as indicated by 
the left hand side of eq. 43 to produce the observed fringe order. This 
is a direct consequence of eq. 37 and the optics of a darkfield polari- 
scope. Knowing the material constants q-j j and the residual effects, eq. 
43 can be used to determine the fringe order produced by a certain 
stress state. It Is important to note that in eq. 43 the five optical 
constants q^- appear as differences in two of the coefficients of 
stress. In reality, only (q^ - q 2 i)» (qi2 " <l22)» and <166 are neces “ 
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sary for determining the effects of applied stress on fringe order. 

This constitutes three quantities Instead of five. In addition, only 
the difference In principal values of the residual dielectric tensor is 
involved, along with the principal residual directions. 

It is legitimate at this point to be concerned with the fact that 
to use eq. 43, such unfamiliar physical constants as the average total 
index of refraction, nl w . the stress-dielectric coefficients, q*i 9 the 
residual dielectric effects, etc. have to be determined. In the next 
section many of these constants will be lumped together into stress- 
optic coefficients and another tensor called the birefringent tensor. 
Then instead of having to be concerned with all of the quantities in eq. 
43, concern will be with, for example, the slopes of the fringe order 
vs. applied stress relations. 


Birefringence Tensor 

To examine other consequences of eq. 37, it is convenient to define 
several quantities which can be thought of as components of another 
tensor. These components combine tensorially to produce the fringe or- 
der Nj. Define 


or 
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where a^- = 
Also define 


2Xn Ive e o 


q-r j i, j - i. 2 , 6, 


'lR 




e lR and N 2R 


T a 2R * 
2xn ave e o 


(47) 


The various quantities N have been called components of the birefringent 
tensor (as opposed to components of the dielectric tensor). Having par- 
ticular values of the components , N^, and is having a particu- 
lar jtate_of_^ij 2 ef^ Using eqs. 46 and 47, eq. 43 becomes 


N T " a 21^°l * ( a 22 " a 12^°2 + ^1R" ^2R^ COS ( 20 R ) } 2 

+ t 2a 66 T 12 + ^1R ‘ ^2R^ sin ( 2 V} 2 

Using eq. 45 this becomes 

N T = " N 2m) + ^1R" ^2r) C0S ( 2e R )} 2 

+ ^^12M + ^1R “ S1 ' n (20p)} 2 


(48) 


(49) 


This equation indicates that the observed fringe count Ny is due in part 
to the effect of normal stresses, N^ M and N^, the shear stress, N^f 
and the residual birefringence. The equation indicates that the fringe 
productions due to normal stresses, N^ M and n£ m , subtract algebraically 
and the difference is added vectorial ly to the fringe production due to 
the shear stress, 2Nj 2M * Residual effect fringe production adds alge- 
braically to the fringe production due to applied load. Equation 49 
states that all effects contribute to the observed fringe. In standard 
polari scopes, the individual components, ne. NJ M , N| m> . ..etc., cannot 
be observed. Only the total effect, Ny, is observable. In a sense. 
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then, the quantities N^, N^.-.-etc. are simply convenient abstractions 
which essentially replace the dielectric tensor. 

Mohr's Circle of Birefringence 

Equation 49 has an interesting physical interpretation. Consider, 
for the moment, the case where the residual effects are zero. Equation 
49 becomes 

N T = ^A N 1M ■ N 2m) 2 + { 2N W 2 ( 5 °) 

If a Mohr's circle is constructed with N^ M and n£ m on the horizontal ax- 
is (i.e., where a and a are located on Mohr's circle of stress) and 
x y 

M? 2 M on the vertical axis, then the state of birefringence due to the 
applied load can be located on these axes. The state of birefringence 
is shown as closed circles on fig. 12a. By geometry, the quantity 

\f [ - 1M 2 S 2 + K 2M ) 2 (51) 

is the radius of the circle. Then twice the quantity in eq. 51 is the 
diameter of the circle. Thus, by eq. 50, with no residual effects pres- 
ent, the diameter of Mohr's circle of birefringence is the fringe order, 
Ny, seen in a polariscope. Equation 50 is the form originally intro- 
duced by Sampson [10] for utilizing the idea of a birefringent tensor. 

He essentially postulated eq. 50 and worked from there. Figure 12b 
shows other important aspects of Mohr's circle of bi refringence. 

A Mohr's circle interpretation can be implied even in the presence 
of residual birefringence. This circle is Mohr's circle of total bire- 
fringence. Referring to eq. 40, define 







( 52 ) 



Construct Mohr's circle by locating nJ m and on the horizontal axis 
and nJ 2H on the vertical axis. The diameter of the circle is 



(53) 


Using the definitions of eqs. 45-47 and substituting into eq. 53, the 
diameter expression becomes 

■ N 2m) + *Ar ' cos ( 26 r )} 2 + 

t^ N 12M + ^1R " si 20^ ^ 

Thus an interpretation of eq, 49 is that the diameter of Mohr's circle 
of total birefringence is equal to the total fringe order, Nj, observed 
in the polariscope. Mohr's circle of total birefringence is illustrated 
in fig. 13. Because of the connection with Mohr's circle, it is obvious 
that the observed total fringe order Ny is also equal to the difference 
in principal values of total bi refringence, i.e. 


where 


N 


1T " oyj * 1T 
2Xn ave e o 


N 1T ■ N 2T’ 

(55) 

and N2T * ? ln T - e 2T" 
2Xn^ 

(56) 

ave o 

(57) 
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Fig. 13 Mohr's circle of total birefringence. 
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When there are no applied loads, N^ M = n| m = = 0- In th ^ s case the 

diameter of Mohr's circle, and the value of the fringe order seen in the 
polarlscope. Is |N 1R - N 2R |. This is the difference in the principal 
values of bi refringence. Also for this no-load case, as will be seen, 
the isocline observed in the polari scope is e R . Thus, the no-load case 
can be used to calculate several important parameters regarding the 
material's optical behavior. 


Redefinition of Coefficients 

If the following definitions are used 
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a ll ' a 21 
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N R 


A A 



(58) a 
(58)b 
( 58) c 
(59) 


eq. 48 becomes 

N T ^ + N r cos( 20 r )} 2 + {S 7 N r sin(2e R )}^ . 


With no residual effects, eq. 60 becomes 




(60) 


( 61 ) 
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This is another form proposed by Sampson [10]. He introduced the three 
stress-optic coefficients fj, f£, and f^* These coefficients can be 
determined for a particular material by applying a known stress and 
observing the fringe order. Specifically, if = 0 , 0 2 = 0, and 
t 12 = then f l can be com P uted - °i ~ 0, o 2 = 0 , and x 12 = 0, then 
f 2 can be determined. If = 0, a 2 = 0, and t^ 2 = x, then f^ can be 
calculated. Alternatively, knowing fj and f 2 and performing an off-axis 
tensile test, f^ can be computed. 

For the case of residual birefringence, calibration of the material 
is more involved. With a no-load condition, N R and 0 R can be deter- 
mined. Then with various known stress states, fj, f 2 , and f 12 can be 
computed. 

For a general case, knowing N R , e R , fj, f 2 » and f^* the fringe 
order can be predicted from a given stress state, 0 ^, o 2 , and t^ 2 , using 
eq. 60. This is exactly what will be done later to predict the fringe 
patterns in the disk. 

In addition, for a given stress state, the optical isocline, 0 T , 
can be predicted. Using eqs. 31e and 35, tan(2e T ) becomes 


tan(20 T ) 



+ 2e 


R 

12M 



(62) 


With the aid of eq. 28a-c this becomes 


tan(20 T ) = 


2e 12M + ( e lR ‘ e 2R^ S, "( 20 R> 


(e 


0 

1M 


- e 


2M 


) + ( e lR - e 2R) cos( 20 r ) 


(63) 


Multiplying the numerator and denominator 
and 47 results in 


by 


2xn ave e o 


and using eqs. 
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tan(20y) = 


2N 12M * ^1R ' ^2r) sin ( 29 R) 
("lM - &) + (N 1R - N 2R ) cos( 26 r ) 


(64) 


Using eqs. 45 leads to 


tan(20y) = 


2a 66 T 12 + ^1R " ^ s1i1 ( 29 r) 


A A 


(<*n “ a 21^1 ” ( a 22 ” a 12 ^ a 2 + t^lR ” ^2R^ cos(20|^) 


(65) 


The definitions of eqs. 58 and 59 produce the final form, namely. 


tan(20y) = 


2 lii + N s1n(2e„) 
'12 

7^ ' + N R cos(2 V 


( 66 ) 


This equation relates the optical isocline to the residual effects, the 
material constants fi, f2* and f^, and the stress state. The optical 
isocline is the parameter observed in the polari scope. It is important 
to note that the optical isocline is a function of the stress level I 
Doubling each stress, for example, changes 0y. The mechanical isocline, 
as © a is often called, can be related to stress from eq. 66 by setting 
residual effects equal to zero. The result is 


tan ( 20^) = 



(67) 


In the presence of residual birefringence, 0^ will never be observed. 
Its existence can only be hypothesized. In the absence of residual ef- 
fects, the mechanical and optical isocline coincide. The mechanical 
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isocline is not a function of stress level. Doubling each stress has no 
effect on tan(20 a ). 

Finally, with no applied stress eq. 66 indicates 

N R sin(26 R ) 

tan < 2 V = N* cos(2e, T = tan < 68 > 

or, as stated earlier, 

8 t = 0 R . (69) 

The observed isocline with no applied load is the residual isocline, 

V 

With the material calibrated for a given stress state, the observed 
isocline can be predicted using eq. 66. This prediction, along with the 
prediction of Ny, will be done for the disk. 

Determining Stresses from Optical Data 

To this point the equations have been manipulated so that if the 
stress state is known, the optical response of the material can be de- 
termined. Specifically, eq. 60 relates the observed fringe order with 
the stress state and residual birefringence. Equation 66 relates the 
observed isocline to the stress state. To determine the stress state 
from the optical data, the equations must be used differently. As with 
isotropic photoelasticity, an auxiliary condition such as the shear dif- 
ference method or a finite-difference form of the two plane-stress 
equilibrium equations is needed to completely determine the stress 
state. However, photoelastic data from polari scope observations pro- 
vides part of the needed information. There are several ways to 


47 


determine just what information is available from a polariscope in the 
way of the state of stress in fiber- reinforced composites. The approach 
taken here starts with eq. 37. 

In eq. 37, subtract the 2nd equation from the 1st and multiply the 
3rd equation by 2. The result is, with a slight rearrangement 

( q ll ' q 21^ a l ‘ ( q 22 ‘ q 12^°2 = ( e lM " e 2H^ ” ( e lM ' e 2H^ 

2q 66 T 12 = 2e 12M ‘ 2e 12M * ( 71 ) 

Using eqs. 28a-c and 31a-c in the right hand sides lead to 


A A A A 

(^11” ^21^1 ” ^22” ^ 12 ^ a 2 = ( e lT“ e 2T^ cos (2®^) “ ( e iR - e 2R^ cos (26p) 

(72) 

2c *66 T 12 = ( e lT ’ e 2T^ sin (20y) ■ “ e 2 R) si n ( 26^) (73) 

It is obvious at this point that multiplying any dielectric tensor 

quantity by £ defines a similar birefringent tensor quantity. 

2Xn ave e o h 

Thus multiplying both sides of eqs. 72 and 73 by T and using eqs. 

2Xn' e 

46, 47, and 55-59 results in 


a-j Op 

J- - J- = N t cos( 20 t ) - N r cos( 26 r ) 


and 


Tip Nj Np 

sin(29 T ) - y - sin(28 R ) . 


(74) 

(75) 


For isotropic photoelasticity, with no residual effects. 
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(76) 


and 

n r - o = e R , (77) 

C being the material’s stress-optic coefficient. Then eqs. 74 and 75 
simplify to 

a - cr v - C N cos(26) (78) 

x y 

t xy = ^ sin(28) , (79) 

where the x-y system has been substituted for the 1-2 material system. 
These are the classical equations of isotropic photoelasticity. 

Equations 74 and 75 indicate that to determine the stress state 
from the optical data, the difference in the principal values of residu 
al birefringence, N R , and the direction of residual birefringence, e R , 
must be known. As stated previously, both these quantities can be 
determined in a polari scope with a no-load condition. 

Transformation of Equations to Global x-y System 

In reality, it is not always convenient to reference the stresses, 
the residual effects, and the material constants to the material system 
Equations 74 and 75 can be transformed to a global x-y system which 
makes an angle of <{> with the 1M-2M system as shown previously. The 
angle <j> is positive as measured counter-clockwise from +x axis to the 
+1M axis. With this definition of <t>, the stresses in the material sys- 
tem are related to the stresses in the x-y system as follows: 
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"cos 2 * si n 2 4> 2s1n$ cos<(» 

$1n 2 $ cos 2 <j) - 2s1 n<|> cos<|> 

•sin* cos$ si n<j> cos<|> (cos 2 <|> - si n 2 <|>) 



• (80) 


Substituting for a^ 9 o^, and in eqs. 74 and 75 results 1° 


cos 2 $ o + si n 2 4> o + 2sin<j> cos<|> t 

* y xy 


$1n 2 $ o v + cos 2 <j> a v - 2sin<j> cos<f> t vw 
- ^ ^ = N x cos( 29 t ) (81) 


and 

^•s1n$ cos$ a + si n<j> cos$ a + (cos 2 <j> - si n 2 <t>)x 

. i ... 

f 12 


- N r s1n(2e R ) 


N N 

•g— s1n(20y) - -j— - sin(20 R ) • 


Gathering terms results in 


( 82 ) 


* 


2(j- + j - ) sin<|> cos* x xy = N t cos( 28 t ) - cos (2e„) , 


(83) 


and 


^i n t co s i^ + ( s1n» cos^ + ( cos**- sin* *' , 
'12 x ' 12 y '12 xy - 

N T . N 

Y~ sin(20y) - -)g~ sin(20j^) • 


(84) 


If the following definitions are made: 
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1 

C 1 = 

j-cos 2 4 » sin 2 <j>j 
f l f 2 

(85) a 

It 

rCos 2 <t> sin 2 <j>^ 

tf 2 ' ^ 1 ' 

(85)b 

1 

2(4 — + -i— ) si n<t> cos<j> 
1 2 

(85)c 

1 _ 

sin* cos<j> 
‘12 

( 85) d 

1 . 

cos 2 <j> - si n 2 4> 
f 12 

(85)e 


then the stress-optic equations become 


+ £-"■ - Ny cos(20y) - N R cos(20 R ) 
12 3 


and 


°x - 3 


v T xv N T N R . 

^ * + = sin(20y) - y ~ sin(20 R ). 


Couple these with the plane stress equilibrium equations 


( 86 ) 

(87) 


and 




8t xv 

_JS£. - 0 
3y 

(88) a 

8a v 

3y~ ~ °* 

(88)b 


and a complete determination of the stresses is possible. Since <j>, fj, 
f£, and are known, - C 5 are known. Note that the quantities ob- 
served in the polari scope, Ny, N R , 0y, and 0 R , enter directly, indepen- 
dent of the x-y system. With certain situations, I/C4 and/or I/C5 might 
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be zero * This Implies C4 = « and C5 = «». If eqs .< 86 and 7 87 are being 
implemented on a computer. It is more convenient to define- • ^ . 


r - 1 
Cl *1 

(89)a 

C 2 = ^ 

(89)b 

C - 1 

C 3 

( 89) c 

r - 1 
C 4-T^ 

( 89 ) d 

r - 1 
C 5 

(89)e 


so that the governing equtions become 

^1 a x ~ ^2 °y + ^3 T xy = c °s(20y) “ cos (26^) (90) 

N t N r 

£4 a x “ ^4 °y + ^5 T xy = 7“ s1n(20y) - ■£— sin(20^) # (91) 

This circumvents having transformed stress-optic properties becoming 
infinite. It is interesting to note that when <|> = 45°, * 0 and the 

equation usually associated with computation of the shear stress, eq. 

91, contains no information about the shear. Also note that for certain 
materials the values of f^, f2» and <|> may combine to produce zero values 
for or C^. 
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In a later section, eqs. 86 and 87 or 90 and 91 will be used to 
determine the stresses In a disk using the photoelastic data and a 
flnlte^dlfference representation of eqs. 88a and b. 
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CALIBRATION OF THE MATERIAL 


Uniaxial Tests 

From the previous section it is clear that to predict Ny and 0y for 
a given stress state, the material properties f 1# f 2 , f^ 2 , Nr and e R 
must be known. Conversely, to predict the stress state from Ny and 0y, 
these same material properties must be known. In the present work, fol- 
lowing the approach of previous investigators [e.g. 11], an uniaxial 
loading of three specimens was initially used to determine fj, f 2 , f^* 
Residual birefringence was not even considered. Thus N R and e R were im- 
plicitly set to zero by using equations which didn't account for residu- 
al effects, for example eqs. 61 and 67. Three calibration specimens 
were cut from the same unidirectional material that was cut into a disk 
for later testing. The calibration specimens were cut as follows: one 

parallel to the fibers (0°), one perpendicular to the fibers (90°), and 
one at 45° to the fiber direction. The material used in these experi- 
ments was manufactured by the Illinois Institute of Technology Research 
Institute under contract with the NASA-Langley Research Center [11]. 

The material had a volume fraction of fibers of 55-60%. This was deter- 
mined by the acid ingestion method with small pieces from the uniaxial 
specimens. The specimens were strain gaged, to determine elastic prop- 
erties, and then subjected to uniaxial loading. The fringe order was 
observed as the load level was varied. For the tensile specimens, at a 
given load, the entire specimen should have been in the same state of 
bi refringence. Bending effects and other eccentricities were minimized 
and felt to be negligible. Yet, during the uniaxial tests it was not 
possible to have the entire specimen in what was felt to be a uniform 
state of birefringence. It was slightly lighter or darker toward the 
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top of the specimen or to one side, never 100% uniform. This was par- 
ticularly true for the higher load levels. Two observers were used to 

minimize val ue- judgement biases in deciding at what load level the spe- 

cimen appeared to exhibit, say, the first fringe. Observations were made 
while loading, while unloading, on different days, and for different 
specimen orientations in the loading frame. At the time the testing was 
being conducted, it was felt local variations in volume fraction or in- 
dex of refraction properties of the fibers and matrix caused the slight 
spatial non-uniformities. This idea was reinforced somewhat by the fact 
that independent of orientation in the loading frame, the same parti ciu- 
lar region of the specimen always exhibited slightly different birefrin- 
gence effects. For the high load levels in the 0° specimen and for all 
load levels in the 90° and 45° specimen, the material was such that the 

fibers were clearly at a different level of birefringence than the 

matrix. 

Figure 14 shows the 0° specimen with no load in the lightfield 
polariscope arrangement. The specimen was 0.625 in. wide, about 12 in. 
long, and 0.090 in. thick. In the photograph the strain gages are 
visible, as are the spots of adhesive securing the lead wires. 

Excessive strain gage adhesive caused the opaqueness in the region 
around the strain gages. There was a long hair-sized flaw along the 
length of the specimen at the bottom end. This is visible in the 
photograph. This flaw did not appear to affect birefringence in the 
region around it. Light centerlines scribed on the specimen are visible 
between the two strain gages. The intersection of these lines was 
chosen as the point at which to count the fringe number. As stated 
earlier, the birefringent effect was not spatially uniform on the 
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Fig. 14 The 0° tensile specimen in light field with no load. 
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specimen and so for consistency the fringe count was always evaluated at 
this same specific point. In the figure the central part of the speci- 
men appears lighter than the upper and lower portions. This is because 

+ •. ’ •'"> v . ‘"'o' v 

the polari scope 1 i ght$ourGe ; and the specimen centerline were aligned on 
the optical axis of the ; polari scope and the light source was a bit too 
intense. This is not related to the aforementioned slight lack of uni- 
formity of birefringence in the specimen. Figure 15 shows the 0° speci- 
men in the light field at what was judged to be a fringe order of 1/2. 
Figure 16 shows the 1st fringe in the light field. Generally, the eye 
discerns spatial variation in light intensity differently than photo- 
graphic film. Because of this, except for the aforementioned difference 
in intensity in the central portion of the specimen, the birefringence 
of the specimens in the photographs presented so far has appeared quite 
uniform and is representati ve of the uniformity until high load levels 
were reached. Figure 17 shows fringe order 4 in the dark-field polari- 
scope. The lack of uniformity of the birefringence is evident in this 
photograph. Bundles of fibers could be clearly distinguished and there 
was a region in the specimen at the lower left which always had a mark- 
edly different level of birefringence. This Is the lighter region in 
the photo. If this had been due to bending effects, there would have 
been an accompanying light region at the top of the test specimen. The 
concern with lack of uniform birefringence was as follows: With a spa- 
tially varying stress state there is a spatial variation of birefrin- 
gence. If the stress st;ate. is spatially uniform and the photoelastic 
sensitivity of the mat^iaT varies spatially, then there is also spatial 
variation in birefringence. If both the state of stress and the photo- 
elastic sensitivity vary spatially in an unknown manner, it would be 
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Fig. 15 The 0° tensile specimen in light field with fringe order 1/2. 
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Fig, 17 The 0° tensile specimen in dark field with fringe order 4. 
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difficult, if not impossible, to extract information regarding the 
stresses from the fringe pattern. Figure 18 shows the 4th order fringe 
in the light field. The lighter region at the lower left of the speci- 
men in fig. 17 appears slightly darker In the situation of fig. 18. 

Figure 19 shows the fringe order as a function of stress for the 0° 
specimen. The experimental data are represented by the solid circles. 
The results of several tests on different days using two observers are 
lumped together. If it Is assumed that there was no Initial birefrin- 
gence in the material, then when - 0, N T is assumed to be zero. For 
the case of no initial birefringence, the hypothesis regarding the ef- 
fect of stress on the dielectric tensor ultimately resulted in eq, 61. 
With a 2 = t 12 =0 in eq. 61, the relation between N and is a simple 
straight line through the origin. Assuming no initial birefringence, 
i.e. including the point = 0, N = 0, a least-square linear fit of the 
data from the 0° specimen produced the solid line in fig. 19. This line 
has a slope of 0.121 x 10~ 3 fringe/psi and an intercept of 0.08 fringe. 
The slope led to a value for f^ of 744 psi/fringe/in. If, as was done 
after seeing these data along with the data from the 90° specimen, it is 
assumed that there was some initial birefringence, then when a j = 0, N-p 
is not necessarily zero. This nonzero value of Np at oj = 0 is really 
N r . The value of Nr could be estimated from the data of fig. 19 by not 
using the points = 0, Ny = 0 In the least-squares analysis. Not in- 
cluding this data point resulted in the linear least-squares fit dotted 
line in fig. 19. The slope of this line is 0.119 x 10~ 3 fringe/psi and 
it has an intercept of 0.13 fringe. The slope led to fj = 756 psi /fringe 
/in. The calibration thus depended on whether or not initial birefrin- 
gence was assumed. It should be pointed out that if at the onset of 
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Fig. 18 The 0° tensile specimen in light field with fringe order 4. 
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Fringe number. 



calibration residual birefringence is suspected, the Tardy method can be 
used to determine the level. 

Figure 20 shows the longitudinal stress-longitudinal strain rela- 
tion for the 0° specimen while fig. 21 shows the lateral strain-longitu- 
dinal strain relation for that specimen. From least-square fits. 

Young's modulus in the fiber direction, Ej # for the specimen was 5.4 x 
10 6 psi and the major Poisson's ratio, was 0.3. 

Testing the 90° specimen was more difficult than testing the 0° 
specimen. The ultimate tensile strength of the material in the matrix 
direction was low compared to the tensile strength in the fiber direc- 
tion. In addition, there was some variation in the ultimate tensile 
strength in the matrix direction from one point to the next in the ma- 
terial. Thus it was not always possible to obtain the same maximum 
fringe count in replicate 90° specimens. A 90° specimen was tested in 
tension until failure and then a portion of the failed specimen was tes- 
ted in compression. Figure 22 shows the set-up used for testing the 90° 
specimen in compression. The loading frame was a hand-operated screw 
device. The 1000 lb. load cell used to measure the compressive load is 
shown. The specimen was strain gaged and it was restrained from 
buckling by sandwiching it between two pieces of plexiglas. The speci- 
men was 0.625 in. wide, 3.6 in. long, and 0.90 in. thick. There were 
holes in the plexiglas restraints so that the restraints were not in 
contact with the strain gages. The plexiglas was not being compressed 
even though the specimen was and so the fringes observed in the specimen 
through the plexiglas were the true fringes in the 90° specimen. 

Figure 23 shows the no-load, lightfield state of birefringence in 
the specimen. Figure 24 shows the 1/2 order, lightfield fringe 
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Fig. 22 Fixture for loading 90° specimen in compression. 
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Fig. 23 The 90° compression specimen in light field with no load. 


(53 


I 


condition. Figure 25 illustrates the darkfield state of birefringence 
at a fringe order of about 1.25. The spatial nonun if ormity of the bire- 
fringence is quite evident here, with nonuni formlty both widthwise and 
lengthwise. As one moves from the top of the specimen toward the bot- 
tom, the image of the birefringence is darker near the centerline of the 
specimen than it is near the edges. This Indicates the center of the 
specimen was at a different level of birefringence than the edges. 

About one-half way down, just above the lower (lateral) strain gage, the 
birefringence suddenly becomes more nearly uniform with width. This 
phenomenon is felt to be related to inhomogeneities in the material 
rather than inhomogeneities in the stress state. In the 90° direction, 
the orthotropic photoelastic material crushes quite easily. Thus any 
deviations from flatness of the ends of the specimen, where the compres- 
sive loads are transmitted to the specimen, would be eliminated by the 
application of load. In addition, in the 90° direction end effects at- 
tenuate quite rapidly due to the low stiffness. Thus any deviations 
from flatness or other end effects would be very unlikely to cause 
abrupt changes in the widthwise uniformity of bi refringence at the cen- 
ter of the specimen. The halos aound the hexagonal bolts were due to 
stresses induced in the plexiglas when the holes were drilled. With the 
situation depicted in fig. 25 it is evident why at times it was diffi- 
cult to determine exactly what fringe order was being observed. For 
this 90° specimen, the fringe order at the center of the specimen was 
used as the fringe order for calibration. 

Figure 26 illustrates the relation between applied stress and 
fringe order for the 90° specimen. Note both tensile and compressive 
stresses are shown. As mentioned previously, Ny is being considered as 
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Fig, 25 The 90° compression specimen in dark field with fringe order 
1.25. 
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a positive quantity and so fig. 26 should actually be plotted as shown 
schematically in fig. 27. However, for purposes of conveniently using 
least-square techniques and for strictly asthetic reasons, plots such as 
in fig. 26 will be used. 

If it is assumed that there was no residual birefringence, then 
when a 2 * 0, the fringe order is zero (see eq. 61 with = T 12 = 0). 

The solid line in fig. 26 shows the least-square calibration with the no 
residual stress assumption. This line has a slope of 0.237 x 10" 3 
fringe/psi and intercept of 0.09 fringe. The slope led to a value of f 2 
= 380 psi/fringe/in. If it is assumed that residual birefringence 
existed, the least-square data fit produced the dotted line. This line 
has a slope of 0.227 x 10~ 3 fringe/psi and an intercept of 0.15 fringe. 
From this slope f 2 = 396 psi/f ringe/inch. Reiterating an earlier dis- 
cussion to avoid ambiguity, the solid line was obtained by using the 
points o 2 = 0, Nj = 0 in the least-square analysis. The dotted line was 
obtained by omitting those data in the least-squares analysis. The 
former approach tends to force the line through the origin when in fact 
it should not go through the origin. 

With both figs. 19 and 26 in hand, residual birefringence was sus- 
pected. With fig. 19 alone it was not clear whether the nonzero inter- 
cept on the N axis was due to the data scatter or due to residual bire- 
fringence. The nonzero intercepts of both figures were, however, strong 
arguments for residual effects. If there were residual effects, were 
figs. 19 and 26 consistent? To determine this, consider eq. 60. 

With <*£ = = 0* eq. 60 becomes 

N T = VM + N r cos( 29 r )} 2 + {N r sin(2e R )} 2 . (92) 
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= 0, eq. 60 becomes 


With Oj = t^2 


(93) 

With an actual polari scope measurement it was quickly established that 
for the uniaxial material, e R was very close to zero. Thus eqs. 92 and 
93 become 



Figures 28a and b illustrate the N vs. a relations for these two equa- 
tions. The presence of residual birefringence, with the residual iso- 
cline aligned with the fibers, results in shifts in the calibration re- 
lations. The 0° calibration is shifted in a direction opposite to the 
shift of the 90° specimen. The experimentally observed shifts in the 
two specimens were consistent with the notion of residual birefringence 
being present, particularly the sense of the shifts. 

Figures 29 and 30 illustrate the elastic properties of the 
90° specimen. Figure 29 indicates that Young's modulus in the matrix 
direction , E 2 , was 1.78 x 10 6 psi and fig. 30 indicates the minor 
Poisson's ration was 0.094. 

A third tensile specimen, with the fibers at 45° relative to the 
tensile direction, was used to complete the calibration. The data from 
this third specimen could be used in several ways. The use depended on 
whether or not residual effects were felt to be important. Figure 31 
shows the no-load state of birefringence in the 45° specimen in the 
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Fig. 28 N t vs. a in the presence of residual birefringence. 
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lateral strain), microinches/inch 




light field. Back-to-back strain gage rosettes were used to measure the 
elastic response of the specimen. Figure 32 shows the lightfield 1/2- 
order state of birefringence of the 45° specimen. This was a relatively 
low load level and yet the fibers caused striations in the birefrin- 
gence.. This could have been due to very high shear stresses at the 
fiber-matrix interfaces. Figure 33 shows the 45° specimen in the dark 
field at what was judged to be the 1st fringe order. As can be seen, it 
was quite difficult to determine what constituted a specific level of 
birefringence. For this specimen the striations were present at all 
load levels. The specimen did, however, go through lighter and darker 
states of 'striated' bi refri ngence. The lightest or darkest states, as 
the load was varied, were used to define the fringe order. 

For the 45° specimen the combined stress state was related to the 

applied stress, a, by the following: 

«1 « o 2 - f ( 96) a 

t 12 = 7 (96)b 

If residual effects were ignored, these stresses could be substituted 
into eq. 61 to yield 

N = a \lj ( 7 — T ~) 2 + • ( 97 ) 

V I 2 12 

The slope of the fringe order vs. applied stress relation, N vs a, can 
be used along with previously obtained values of f^ and f 2 to determine 
the value of f 12 . Figure 34 shows the fringe order-appl ied stress rela- 
tion for the 45° specimen and also indicates a least-square straight 
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Fig. 32 The 45° tensile specimen in light field with fringe order 1/2. 
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line. Since it was assumed in fig. 34 that residual effects were not 
important, the points a = 0, N T = 0 were used as data in the least- 
square scheme. The slope of the line, 0.288 x 10 3 fringe/psi, resulted 
in a value of 312 psi/fringe/in for f^. 

If it is assumed residual effects were important in the 45° speci- 
men, eq. 60 must be used. It is more convenient to square eq. 60 to 
obtain 


2 

n t 

-[} (^-f-) 2 

4 T i 

+ (t -) 2 ]* 2 + [Mr- 
r 12 K T 1 

• + n r * 

(98) 

where 

eq. 96 has been 

used along with the fact that e R = 0. 

If A, B 

and C 

are defined as 







A = [j (f - - f-) 2 
* T 1 T 2 

* r- 

r-) 2 ] . 
'l2 

( 99 ) a 



6 - V 77 - 


* 

(99)b 



c- n 2 r 

> 


( 99) c 


then eq. 98 is of the form 


Nj * Ao 2 + Ba + C . (100) a 

2 

If the experimental measurements of Nj and a are plotted as Ny vs. a, 
and a least-squares parabola of the form 

Ny = ao 2 + go + y (100)b 

is fit to the data, then A, B, and C can be estimated from 
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A - a (101)a 
B - g (101)b 
C * y . (101)c 


If and f2 have been determined from previous calculations, then a can 

2 

be used to determine f 12 . The value of y should give a value of N R con- 
sistent with the value obtained in 0° and 90° tests. The value of B 

should be consistent with the values of N R , fj, and f 2 . Figure 35 shows 

2 

the relation between Hj and a. The least-squares parabola has the fol- 
lowing coefficients, ~ '•* 

a = 0.75 x 10“ 7 ^f|inge/psi ) 2 (102)a 

3 = 0.45 x 10” 4 (Tringe) 2 /psi (102)b 

Y = 0.049 (fringe) 2 ( 102) c 

This value of a led to « value for f^ 2 of 335 psi /fringe/in. The value 
of y predicted a value of N R = 0.22. This value of N R was higher than 
the value of N R compute|[‘|From the 0° or 90° specimens. However, there 

y 

was much more difficulty in determining the level of bi ref ringence in 

-3 * 

vvl 

the 45° specimen and so, as seen in figs. 34 and 35, there was much 
scatter to the Ny - a data. Hence, stress-optic properties obtained 
from this specimen may not be as reliable as properties taken from the 
0° and 90° specimens. Unfortunately, with this approach to calibration, 
fj 2 musing determined from this 45° specimen. 

In addition to the value of N R being different than previous calcu- 
lations, the value of 3, eq. 102b, was inconsistent. The quantity 3 


85 




approximates B in eq. 100 (see eq. 99b). Since > f 2 , B will be a 
negative number. The value of 3 from the least-squares approximation 
was positive. The reason for this inconsistency was not evident. It is 
not felt to be due entirely to the data scatter responsible for the 
inconsistent value of N R . 

Figure 36 shows a plot of the shear stress t 12 as a function of the 
sheer strain The shear strain was determined from the rosette ga- 

ges while the shear stress was computed using eq. 96b. For some reason 
the strain gage readings were erratic and so the elasticity data had 
scatter. A least-square straight line through the t - y data yields a 
value of G 12 570,000 psi. The longitudinal stress-longitudinal 

strain ralation for this 45° specimen resulted in a Young's modulus, 
i.e. of 1.56 x 10 6 psi. 

Table 1 summarizes the optical properties and Table 2 summarizes 
the elastic properties of the material as obtained from the uniaxial 
tests. 

Four Point Bend Tests 

After it was hypothesized that residual birefringence was present 
in the material, it was decided to use four-point bending tests to op- 
tically calibrate the material. There were several reasons for this. 
First, it was felt the stress gradients in a bending specimen would re- 
sult in sharp, more well defined fringes. It would thus be easier to 
judge the presence of, say, the first fringe order. Using the flexure 
formula and determining the distance from the neutral bending axis to a 
particular fringe, the stress levels producing that fringe order could 
be accurately determined. Second, with residual birefringence in the 
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Table 1 


Optical properties of unidirectional material 
as determined from uniaxial tests 


Specimen 

No Initial 

Initial 

Type 

Birefringence Assumed 

Birefringence Assumed^ 1 ) 

0° 

fj = 744 psl/frlnge/ln 

fj = 756 psi/fringe/in 
N r = 0.13 

90° 

f2 = 380 psl/frlnge/ln 

f 2 = 396 psl/frlnge/ln 
N r = 0.15 

45° 

fl 2 - 312 psi/fringe/in 

f 12 - 335 psi/fringe/in 
N r = 0.22( 2 ) 


(1) 9 r measured to be zero 

(2) N R and e not consistent with previous computations 
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Table 2 


Elastic properties of unidirectional material 
as determined from uniaxial tests 

Specimen Type Elastic Properties 

0° Ej = 5.4 x 10 6 psi v 12 * 0.3 

90° E 2 = 1.78 x 10 6 psi v 21 = 0.09 

45 ° Gl 2 = 0.57 x 10 6 psi E 45 = 1.5 x 10 6 psi 
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material, the fringes would not be symmetrically disposed on either side 
of the bending neutral axis. The 'optical neutral axis' would be off- 
set from the conventional neutral axis, i.e. the centerline of the bend- 
ing specimen. This off-set could be used to determine the residual 
fringe order Nr. Finally, one test could be used for both tensile and 
compressive stresses. With a known applied bending moment and using 
Tardy compensation to determine fractional fringe orders, considerable 
stress-optic data could be obtained from the four-point bending arrange- 
ment. 

All the previously discussed uniaxial specimens were used as four- 
point bending specimens. These specimens were reused as bending speci- 
mens for two reasons. First, it conserved valuable photoelastic materi- 
al. Second, it insured that variability of optical properties from lo- 
cation to location on the panel would not contribute to any differences 
in the calibration constants between the uniaxial tests and the four- 
point bend tests. There was one disadvantage in using the uniaxial 
specimen. These particular specimens, as stated before, were approxim- 
ately 0.625 in. wide. Wider specimens would have resulted in more ac- 
curacy in resolving distances from the neutral bending axis. Actually, 
photographs of the four-point bend test fringe patterns were enlarged to 
aid in determining distances from the neutral axis. Thus the size of 
the specimen was not important from that standpoint. However, using a 
wider specimen as the basis for scaling would have resulted in more res- 
olution than using a narrower one. 

Figures 37-39 show the three four-point bending tests. Figure 37 
shows the 0° specimen. In this case the fibers were parallel to the 
neutral bending axis (horizontal in the figure). Figures 38 and 39, 
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Fig. 38 Four-point bend test of 90° specimen. 





Fig. 39 Four-point bend test of 45° specimen. 


I 


respectively, show the 90° and 45° specimens. With the 90° orientation, 
the fibers were perpendicular to the neutral axis. The 90° specimen was 
somewhat shorter than the other two since it was the remains of a failed 
tensile specimen. Recall, the 90° tensile specimens were tested in ten- 
sion to failure and then one of the remaining pieces was tested in com- 
pression with the buckling restraints. 

With the specimens as shown in figs. 37-39, the portion of the beam 
above the neutral axis was in compression. For the 0° specimen. 





(103) 


where z was measured positive downward from the neutral axis and M and I 
have their usual meaning. Equation 94 provides the relation between 
and Nj, i.e. 



(94) 


Using various values of moment, looking at several z locations, and 
using Tardy compensation, sufficient data could be generated to reliably 
determine fj and N R . Figure 40 shows the relation between and Ny. 

The maximum stress levels used in the four-point tests were not neces- 
sarily identical to the stress levels in the uniaxial tests because in- 
terest here was in residual birefringence. The slope of the line in 
fig. 40 is 0.106 x 10“ 3 fringe/psi and the intercept, N R , is 0.11. The 
slope led to a value of f^ of 849 psi/fringe/ inch. This is 11% higher 
than the value of the fj using the tensile tests, assuming an initial 
bi refringence. 
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For the 90° specimen. 


_ Mz 

°2 = T~ 


(104) 


and eq. 95 provides the relationship between and Nj, i.e. 



Figure 41 shows the relation between Ny and a^. The slope of the rela- 
tionship is 0.234 x 10" 3 fringe/psi and the intercept is -0.07 fringe. 
The slope resulted in a value of f 2 of 385 psi/fringe/inch. This is 
about 3% lower than the value obtained in uniaxial tests. 

For the 45° specimen, eqs. 96 and 98 are applicable, where in those 
equations 


a 



(105) 


Figure 42 shows the relation between and a. A least-squares parabola 
of the form eq. 100b was fit through the data. The coefficients were 


a = 0.837 x 10' 7 (fringes) 2 /(psi ) 2 (106)a 

£ = 0.34 x 10" 4 (f ri nges) 2 /psi (106)b 

Y = 0.0451 fringes . (106)c 

Using the values of fj and f 2 from the 0° and 90° bending tests, the 
value of a led to a value of f 12 of 319 psi/fringe/in. This is about 
10% lower than the value of f^ 2 obtained from the uniaxial tests. The 
value of y led to a residual fringe value of 0.21. The sign of 3 was 
again inconsistent with the other data. Table 3 summarizes the optical 
calibration from the four-point bending tests. 
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tic relation for 90° specimen from four-point bend test. 
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Fig. 42 Stress-optic relation for 45° specimen from four-point bend test. 



Table 3 


Optical properties of unidirectional material 
as determined from four-point bend tests. 


Specimen Type 


90 c 


45 < 


Data( *) 


R = 


12 


849 psi/fringe/in 

0.11 

385 psi/fringe/in 
0.07 

319 psi/fringe/in 

0.2l( 1 2 ) 


(1) e R measured to be zero 

(2) 3 not consistent with previous computations 
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In summary, it Is felt that the four-point bending tests yielded a 
more accurate calibration. The ambiguity concerning fringe definition 
in the uniaxial tests did not occur in the bending tests. Furthermore, 
having tensile stresses and compressive stresses at the same time was 
felt to be valuable. Since the transition between tension and compres- 
sion is associated with determining the value of Nr, having a single 
testing fixture and simultaneously determining the value of N T associa- 
ted with tension and the value associated with compression should reduce 
the variance in the determination of N R . 
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FORWARD TESTING: GLOBAL BEHAVIOR 


With the material calibrated, eqs. 60 and 66 could be used to pre- 
dict the isochromatic and isoclinic fringes in a disk, assuming the 
stresses were known. Alternatively, using eqs. 74 and 75, in addition 
to a third condition, the stresses could be computed, assuming the iso- 
chromatic and isoclinic fringe patterns were known. The former has been 
referred to here as forward testing while the latter has been referred 
to as backward testing. The forward testing is discussed in this and in 
the next section 

The comparing of the predicted fringe patterns with the experiment- 
ally observed fringe patterns was done using two approaches. Both com- 
parisons were done using an orthotropic disk in diametral compression. 
The first comparison was done by looking at the disk as-a-whole. This 
is discussed in this section and is referred to as global behavior. The 
second approach is described in the next section and is referred to as 
local behavior. To study the global behavior, photographs of fringe 
patterns in the disk were compared with computer-generated fringe pat- 
terns. These computer-generated fringe patterns were based on eq. 60 
and the elasticity solution of Appendix B. This solution is based on 
the complex variable approach of Lekhnitskii [12]. The stress-optic 
data required in eq. 60 was taken from the tensile and four-point bend- 
ing tests and the stresses, aj, » and were determined from the 
elasticity solution. Figure 43 shows the set-up used to load the disks 
in diametral compression. The disk used was 3.12 in. in diameter and 
0.090 in. in thickness. Shown in the photograph is the plexiglas sand- 
wiching plates used to prevent the disk from buckling under the 
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Fig. 43 Fixture for loading disks in diametral compression. 
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compressive load. This loading fixture is the fixture described in 

conjunction with the discussion of the 90° compression calibration 

specimens. 

One dilemma in the comparison of theory and experiment was the 
choice of the numerical values of the stress-optic parameters, particu- 
larly N R . The calibration tests produced several values of N R , ranging 
from 0.07 to 0.22. Realistically, the 0.22 value from Table 1 and the 
0.21 value from Table 3 were unreliable. As previously indicated, they 
were not consistent with the other calibration data from the same speci- 
men. It was not clear a simple average of the various values would be a 

valid choice for the value of N R . A value of 0.15 was chosen, somewhat 

arbitrarily, for N R . This certainly was a good approximation to the 
values measured from the tests. The values of the other stress-optic 
coefficients were as follows: 

fl - 849 psi /fringe/in. 
f 2 = 385 psi/fringe/in. 
f 12 = 319 psi/fringe in. 

These data are from Table 3. 

Figure 44 shows a typical computer-generated fringe pattern for the 
disk. For the case shown, the fringe order is Nj = 0.5 and the fibers 
are 30° relative to the loaded diameter. Figures 45, 48, 51, 54, and 57 
show the comparisons between the observed fringe patterns in the disk 
and the computer-generated fringe patterns. The figures show the fringe 
patterns for the disk with the load at various angles relative to the 
fibers. In all figures the primary fringe order is Ny = 0.5. Other 


(107) a 
( 107 ) b 
( 107 ) c 
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Fig. 44 Computer-generated fringe pattern for disk in diametral 
compression, Nj=0.5. 


105 





fringes are visible near the application of the load. These fringes 

were not studied because of the lack of ideal loading and lack of accur- 

\ 

acy in the elasticity solution in this region. Visible in the photo- 
graphs of the disk is crazing damage on the disk edge appearing as dark 
spots on the edge. When the load was applied, the disk material in the 
immediate vicinity of the load application point would craze and thus 
become opaque to light transmission. At first, the load was introduced 
into the disk with a rounded contact point. The radius of curvature of 
the rounded contact was 0.125 in. This small radius was used to simu- 
late point loads. However, the damage induced in the disk was severe. 
Thus the flat contacts seen in the figure were used. Still some crazing 
would always occur, especially when loading the disk perpendicular to 
the fibers. 

Figure 45 shows the observed and predicted fringe pattern when the 
load was applied along the fibers. The load level was 500 lb. Certain- 
ly the characteristics, such as the vertical bow-tie shape, correlated 
well. On the one hand, it is quite remarkable how good the correlation 
was. The state of stress in the disk was quite complex. The stress- 
optic law is a hypothesis and the calibration of the stress-optic coef- 
ficients was done with simple tests. The three main ingredients seemed 
to combine well to lend credibility to the whole notion of orthotropic 
photoelasticity. On the other hand, a closer scrutiny of the results 
shows the predicted horizontal distance from the center of the disk to 
the Nj = 0.5 fringe was less than observed. Scaled to the disk's actual 
size, the distance from the center to the fringe was predicted to be 
0.33 in. From the photograph the distance was closer to 0.40 Jn. In- 
terestingly enough, if it was assumed the residual fringe value, N R , was 


106 




Fig. 45 Comparison between predicted and observed fringe pattern, 
fibers aligned with load, N T =0.5 (predicted=white dots). 
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0.07, one-half the value being assumed, the distance prediction was very 
close to the value observed. Figure 46 shows the predicted fringe pat- 
tern as a function of three levels of residual bi ref ringence, N R = 0, N R 
= 0.07, and Nr - 0.15. The level of residual birefringence obviously 
makes a difference in the fringe pattern generated in the disk. Figure 
47 shows another interesting aspect of the effects of residual bire- 
fringence. The figure shows the predicted effects of simply reversing 
the load and subjecting the disk to diametral tension. For no residual 
birefringence, the sign of the load has no influence on the observed 
fringe pattern. With residual bi refringence, the sign of the load 
influences the fringe pattern produced. 

Figure 48 shows a comparison between predicted and observed fringe 
pattern in the disk when the load was oriented 30° relative to the fi- 
bers. The load level was 325 lb. (With the sign convention adapted in 
Appendix B, this corresponds to a = -30°. ) The fibers caused the 
streaking in the photo and also contributed to the lack of sharp fringe 
definition in some regions of the disk. This obscuring was the same as 
observed in the testing of the uniaxial specimens described earlier. 

The distance from the disk center to the fringe location, along the di- 
ameter perpendicular to the fibers, was slightly overpredicted by the 
theory. Figure 49 shows the effect of various levels of residual bire- 
fringence on the predicted fringe patterns. The differences are subtle 
but there are slight differences. When shear stresses, in the 1M - 2M 
coordinate system, dominate the stresses, small changes in residual bi- 
refringent effects become less important. This is because for this ma- 
terial e R = 0 and only the term involving and o? in eq. 60 is affect- 
ed by varying Nr. Figure 50 shows the effect of a load reversal on the 
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Fig, 46 Effect of level of residual bi refringence on predicted fringe 
pattern, fibers aligned with load, Nj=0.5. 
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Fig. 47 Effect of residual birefringence on the load reversal fringe 
pattern, fibers aligned with load, N-pO.5. 
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Fig. 48 Comparison between predicted and observed fringe pattern, 
fibers 30° relative to load, Ny=0.5 (predicted=white dots). 
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predicted fringe pattern for this 30° orientation. Again, it is clear 
that simply reversing the load direction has an impact on the fringe 
pattern. 

Figure 51 shows the predicted and observed fringe patterns for the 
case of the fibers making a 45° angle with the diametral compressive 
load. The load level was 325 lb. (This orientation corresponds to a = 
+45°.) Again the fibers tended to spread the fringe but the pattern was 
quite close to the predicted configuration. A close scrutiny of the 
distance from the disk's center to the half-order fringe, along a diame- 
ter perpendicular to the fibers, showed the predicted value to be quite 
close to the observed value. Figures 52 and 53 show the effects of 
residual level and load reversal , respectively. 

Figures 54-56 shows similar comparisons for the case with the 
fibers 60° relative to the load direction. Again the load level was 325 
lb. (This corresponds to +60 in the nomenclature of Appendix B). 

The spreading of the fringe by the fibers was quite severe. However, 
the computer-generated fringe pattern shows good correlation between the 
dots of the predicted fringe pattern and the center of the fringe. 

Figure 57 shows the observed and predicted fringe pattern for the 
situation when the load was 90° relative to the fibers. Here the load 
level was 520 lb. The correlation for this case was the poorest of all 
cases. Using a residual birefringence of 0.07 made the correlation 
better. Figure 58 shows that the predicted fringe pattern for this ori- 
entation of load and fibers is fairly sensitive to the level of residual 
birefringence. As expected. Fig. 59 shows that with residual bire- 
fringence the sign of the load is important in determining the nature of 
the observed fringe pattern. 
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Fig. 50 Effect of residual birefringence on the load reversal fringe 
pattern, fibers 30° relative to load, Ny=0.5. 
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Fig. 51 Comparison between predicted and observed fringe pattern, 

fibers 45° relative to load, Ny=0.5 (predicted=white dots). 
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Effect of residual bi refringence on the load reversal fringe 
pattern, fibers 45° relative to load, Nt=0.5. 



Fig. 54 Comparison between predicted and observed fringe pattern, 

fibers 60° relative to load, N T =0.50 (predicted=white dots). 
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Fig. 57 Comparison between predicted and observed fringe pattern * 

fibers perpendicular to load, N-pO.5 (pred1cted=wh1te dots). 
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Fig. 59 Effect Of residual birefringence on the load reversal 
pattern, fibers perpendicular to load, Ny=0.50. 


| R = °- 15 
I -P, n r = 0 

N r = 0.15 


fringe 
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In conclusion, it appears that the idea of orthotropic photoelasti- 
city is, on-the-whole, a valid concept. However, fringe definition can 
definitely be a problem. In the cases here, the ‘answer* was known and 
so it is easy to conclude the fringes were, more or less, where they 
should have been. However, for the 60° orientation (fig. 54), at the 
two o'clock and eight o'clock locations on the disk (high noon being at 
the upper load application point), it may not have been possible to 
sketch in the fringes if it were not known a priori where they were sup- 
posed to be located. This is a serious concern in the practical appli- 
cation of orthotropic photoelasticity. This prompted the series of 
tests described in the following section, namely, a closer look at the 
predicted and observed fringes. 
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FORWARD TESTING: LOCAL BEHAVIOR 


Whereas forward testing of the global behavior indicated an overall 
comparison between theory and experiment, a more detailed comparison was 
of value. For this reason another series of comparisons between theo- 
retical and experimental fringe data was made. Specifically, the total 
isochromatic fringe order, Nj, and the optical isocline, Qj, were com- 
pared point-by-point across a diameter perpendicular to the load and 
across a diameter perpendicular to the fibers. These comparisons were 
made for diametral loads at various angles relative to the fibers. 

These comparisons are shown in figs. 60-75. The figures show the exper- 
imentally measured values and their theoretical predictions. On each 
figure there are several theoretical predictions, each prediction corre- 
sponding to a different level of assumed residual birefringence, N R . 

The three levels are the same as in the previous section. Nr = 0, 0.07, 
and 0.15. The isochromatic fringe predictions were based on eq. 60, the 
stress-optic data for the material, and the complex variable elasticity 
solution. The optical isocline predictions were based on eq. 66, the 
stress-optic data for the material, and the elasticity solution. On 
each figure the horizontal axis shows the nondimensional distance across 
a diameter as measured from the left edge of the specimen. Measurements 
were only made from one edge of the disk to the center. The vertical 
axes show the isochromatic or isoclinic measure, depending on the fig- 
ure. On the figures, the open squares represent the experimentally mea- 
sured values. The solid circles represent the theoretical predictions 
for the case of no residual birefringence (N R = 0). The solid triangles 
represent the predictions for the case of Nr = 0.07 and the solid 
squares represent the theoretical values for the case N R = 0.15. 
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Figure 60 shows the comparison of isochromatic fringe order for the 
case when the loading was aligned with the fibers. The load level used 
was 407 lb. Obviously the diameter perpendicular to the load and the 
diameter perpendicular to the fibers were the same for this case. As 
can be seen, starting at the left edge of the disk and moving toward the 
center, the experimentally measured isochromatic fringe count decreased 
to practically zero and then increased to a value slightly greater than 
one-half at the center. Independent of the assumed level of initial bi- 
refringence, the theoretically predicted fringe count also approached 
zero somewhere along the radius. The initial level of bi refringence 
only affected where it went to zero. For no initial birefringence (sol- 
id circles) the fringe count dropped to zero at x/D = 0.15. Experimen- 
tal measurements (open squares) showed minimi a! fringe value at x/D = 
0.21. A residual level of birefringence of 0.07 caused the predicted 
fringe level to drop to zero at x/D = 0.23. The zero fringe was pre- 
dicted to occur at x/D - 0.3 when the residual birefringence was 0.15. 
Except for near the left edge, the experimental observation correlated 
well with the predictions based on low (N R between 0 and 0.07), as op- 
posed to high (N R = 0.15), levels of residual birefringence. 

Figure 61 shows the predicted and observed optical isocline for the 
case of the load aligned with the fibers. This figure is quite inter- 
esting. In reality all symbols should be plotted at 0 = 0° or 0 = -90°. 
To illustrate an important point, the symbols have been shifted slightly 
off the 0° and -90° values. It was the rapid jump of the isocline from 
0° to -90° as one moved in along a radius that was important. The x/D 
location at which the optical isocline jumped from 0° to -90° was the 
same location at which the total isochromatic fringe count dropped 
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Fig. 60 Comparison between predicted and observed isochromatic 
fringe order along diameter perpendicular to load, 
fibers aligned with load. 
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isoclinic angle. 



Fig. 61 Comparison between predicted and observed isoclinic 
angle along diameter perpendicular to load, fibers 
aligned with load. 


to zero. This location can be thought of as an optical isotropic point. 
At such a point the total isochromatic fringe order is zero and the 
isoclinic is indeterminant. Referring to eq. 60, the second squared 
bracketed term is zero along the diameter. There is no shear stress, 
and the residual isocline, 6 R , is zero. With a proper combination 
of normal stresses, a ^ and a^, the first squared bracketed term can be 
made zero. This was precisely the case at the various x/D locations. 

As can be seen from eq. 60, the value of N R affects this term going to 
zero. Examining eq. 66, if the first and second terms of eq. 60 are 
zero, eq. 66 becomes 


tan 20-j- = -§■ • 


(108) 


This leads to i ndetermi nent values for Qj. It should be noted, not only 
does the value pf Nr affect the location of the isotropic point, but the 
load level does also. Since N R is fixed and and can be varied by 
varying the load level, the physical location of the optical isotropic 
point can be controlled. This was not explored in this study. 

From the above discussion, it seems that it may be possible to de- 
termine N^ by comparing the theoretical and experimental isotropic point 
location along the diameter. Since the isotropic point in the experi- 
ment was located between the predicted isotropic point for no initial 
bi refringence and the case N R = 0.07, the assumed lower levels of resid- 
ual birefringence seem to correlate better with experiment than does the 
assumed higher level. 

Figures 62 and 63 show further isochromatic and isoclinic data for 
the disk with the load aligned with the fibers. Since the diameter 
perpendicular to the load and the diameter perpendicular to the fiber 
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Fig. 62 Comparison between predicted and observed isochromatic 
fringe order along 30° diameter, fibers aligned with 
load. 
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Fig. 63 Comparison between predicted and observed isocline 
angle along 30° diameter, fibers aligned with load 




were the same, data was studied along another diameter. This diameter 
was arbitrarily picked to be a diameter which made a 30° angle with the 
horizontal diameter. 

Moving from the edge of the disk toward the center along the 30° 
diameter produced trends in the data similar to the horizontal diameter 
data. The isochromatic fringe number decreased and then increased to 
about 0.5. The minimum value of the theoretically predicted isochromat- 
ic varied, depending on the assumed value of initial bi refri ngence. For 
no assumed initial bi refri ngence (solid circles), the isochromatic 
fringe value dropped to zero at x/D = 0.12. The experimental ly measured 
isochromatic never dropped to zero. The same was true for the theoreti- 
cal predictions based on assumed initial birefringence values of 0.07 
and 0.15. 

Figure 63 shows the theoretical ly predicted and the experimental ly 
measured optical isocline parameter for the 30° diameter. There were 
vast differences, depending on the assumed level of initial birefring- 
ence, between the predicted and the measured values. The experimentally 
measured isocline varied smoothly from about 5° at the outer edge, 
through 0° at x/D = 0.06, to -85° at the center of the disk. The theor- 
etical cases of N R = 0.07 and N R = 0.15 also varied smoothly but the 
case N r = 0 exhibited a sharp jump at x/D = 0.12. This was the same lo- 
cation that Nj = 0 for this theoretical case. This behavior was similar 
to the jump behavior in the optical isocline along the horizontal diame- 
ter for this fiber/load orientation. Again, except for the outer edge, 
the theoretical predictions for the case N R = 0.07 seemed to correlate 
the best with experimental observations. The isochromatic predictions 
almost coincided and the isocline predictions were close. 
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Recall, there was concern over choosing the value of N R from the 
tensile and compression test data. There were values, though unrelia- 
ble, of N R as high as 0.21. There was strong evidence from the tensile 
and compression data that N R was in the range 0.07 - 0.15. The value of 
0.15 was chosen from the test data of the 0°, 90°, and 45° specimens. 

Yet the experimentally observed global and local behavior correlated 
better with theoretical predictions obtained by using a lower value of 
assumed residual birefringence, specifically N R = 0.0 or 0.07. This 
trend toward good correlation with 0 < N R < 0.07 in the global and local 
tests contradicted the findings of the 0°, 90°, 45° tests. The contra- 
diction is still a puzzle. 

Figures 64 and 65 show the isochromatic and the isoclinic behavior 
along a horizontal diameter for the case of the fibers rotated 30° from 
the applied diametral load. The observed isochromatic and isoclinic be- 
havior was similar to the other two cases discussed. The experimentally 
observed isocline changed from negative values to positive values at x/D 
of about 0.08, The theoretical prediction for the case N R = 0 showed a 
very sharp jump in isocline at this x/0 location. At the same time, the 
isochromatic prediction for this level of bi ref ri ngence practically 
dropped to zero at this same location. Throughout the testing of all 
cases presented in the figures in this section, when the isocline exper- 
ienced a rapid change with spatial location, the isochromatic fringe 
number was going to a minimum value. This was true with the experiment- 
al observations and all theoretical predictions. It seemed that the 
closer the isochromatic fringe approached zero, the more dramatic was 
the jump in the isoclinic. 
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Fig. 64 Comparison between predicted and observed isochromatic 
fringe order along diameter perpendicular to load, 
fibers 30° relative to load. 
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Fig. 65 Comparison between predicted and observed isoclinic 
angle along diameter perpendicular to load, fibers 
30° relative to load. 
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Several additional comments should be made regarding the data in 
figs. 64 and 65. In the range 0.2 < x/D < 0.5, the assumed level of re- 
sidual birefringence did not affect the predicted i sochromatic. All 
predictions coincided. Because of this, in that same range along the 
diameter, the correlation between predicted and observed isochromatics 
didn't seem to favor any particular level of N R . However, the isoclinic 
figure showed reasonable correlation with the N R = 0.07 case, particu- 
larly the characteristics near the change in sign in the isoclinic. 

Figures 66 and 67 show more i sochromatic and isocline data for the 
case with the fibers oriented 30° relative to the load. These figures 
are for the diameter perpendi cular to the fiber direction. Observations 
and comments similar to the discussion of previous figures can be made 
regarding these two figures. As seen in fig. 67, for example, the ex- 
perimentally observed isoclinic experienced a sharp change at 
x/D » 0.14. It was at this location N T was a minimum. At first, the 
sharp changes in isoclines were viewed as an undesirable feature. Even- 
tually, however, these sharp changes were a welcome characteristic. 

They served as sort of a datum for minimum i sochromatic fringe counts 
and a point where there was most likely a change in sign in the isoclin- 
ic. Polariscopes, by their nature, respond similarly to isoclinics of 
complementary values, i.e., isoclinics that differ in value by 90°. 

When the sign of the isoclinic changes, or the isoclinic otherwise goes 
to zero, interpretation of the data can be confusing. The existence of 
the sharp changes helped clarify the interpretation. 

Figures 68-75 show the isocline and i sochromatic data for cases 
when the fibers were at 45° and 90 to the load direction. These cases 
will not be discussed in detail. In light of previous discussions, some 
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Fig. 66 Comparison between predicted and observed isochromatic 
fringe order along diameter perpendicular to fibers, 
fibers 30° relative to load. 
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Fig. 67 Comparison between predicted and observed isoclinic 
angle along diameter perpendicular to fibers, fibers 
30° relative to load. 
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Fig. 68 Comparison between predicted and observed isochromatic fringe 
order along diameter perpendicular to load, fibers 45° relative 
to load. 


139 



isoclinic angle 




# theory, N = 0 
R 

A theory, N R = 0.07 
■ theory, N R = 0.15 
□ experiment 


w • • i t 

* * * * * » sc 


i i i 2! I 



Fig. 69 Comparison between predicted and observed isocline angle along 
diameter perpendicular to applied loads, fibers 45° relative to 
load. 
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Fig. 70 Comparison between predicted and observed isochromatic 
fringe order along diameter perpendicular to fibers, 
fibers 45° relative to load. 
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Fig. 71 Comparison between predicted and observed isoclinic 
angle along diameter perpendicular to fibers, fibers 
45° relative to load. 
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Fig. 72 Comparison between predicted and observed isochroma tic 
fringe order along diameter parallel to fibers, fibers 
perpendicular to load. 
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Fig. 73 Comparison between predicted and observed isoclinic 
angle along diameter parallel to fibers, fibers 
perpendicular to load. 
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Fig. 74 Comparison between predicted and observed isochromatic 
fringe order along 30° diameter, fibers perpendicular 
to load. 
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75 Comparison between predicted and observed isoclinic 
angle along 30° diameter, fibers perpendicular to load 




of the observations and comments are obvious. The load level used for 
the 45° and 90° orientations was 350 lb. It can be seen from the 
figures that the experiments correlated better with the lower levels of 
assumed residual bi refringence than they did with the higher levels. 

It was of interest to put these experiments of local behavior into 
context by repeating the experiments with a standard isotropic disk. A 
great effort is made in the manufacture of isotropic photoelastic mater- 
ial to eliminate residual bi refringence. Thus the isotropic tests could 
not be used to check the effects of residual birefringence. However, 
the tests could certainly be used to check the spatial distribution of 
isochromatic and isoclinic data. If it is assumed the theoretical ly 
predicted data were correct, then the experiments would evaluate how 
well photoelastic data was being measured. If it was assumed that the 
photoelastic measurements were correct, then the comparison between ex- 
periment and theory could be viewed as an evaluation of the theory of 
isotropic photoelasticity. The theory of isotropic photoelasticity has 
been verified time and time again and so experiments with the isotropic 
disk should be looked upon as an indication of the accuracy to which 
isochromatic and isoclinic fringe data could be recorded with the 
pol ari scope. 

Figures 76-79 show the comparison between the measured and the pre- 
dicted isochromatic and isoclinic fringe information for the isotropic 
disk. The predictions are based on the exact elasticity solution for an 
isotropic disk as presented in Frocht [13], The load level used was 125 
lb. The isotropic disk material was PSM-1. Figure 76 shows the 
isochromatic fringe numbers across a diameter perpendi cul ar to the load. 
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Fig. 76 Comparison between predicted and observed isochromatic 
fringe order along diameter perpendicular to load, 
isotropic disk. 



Figure 77 shows the isoclinic data for that diameter. The correlation 
was very good. Figures 78 and 79 show similar information for a diame- 
ter oriented at 30° to the horizontal diameter. Again the correlation 
was excellent. From these four figures, there was no question as to the 
ability to make accurate measurements with the polariscope. 

From these studies with the orthotropic and isotropic disks, one 
conclusion is certain. Orthotropic photoelasticty , as it is hypothe- 
sised here, is not as accurate as isotropic photoelasticity. Although 
the conclusion is straight forward, it is not so clear why there is a 
difference in accuracy. It can, however, be traced to two sources. 
Possibly, the hypothesis regarding the stress-optic behavior, given by 
eq. 36, is not entirely valid. The averaging of the effects of the fi- 
bers and matrix, in a through-the-thi ckness sense, may not be as accur- 
ate for optical behavior as it is for, say, elastic behavior. In addi- 
tion, there may be an interaction between shear and elongation effects. 
That is, maybe eq. 36 should not be orthotropic in any coordinate sys- 
tem. These issues are directed at the theory itself. Perhaps the inac- 
curacy comes in the practicing of the theory. Fringe definition in or- 
thotropic photoelastic material is not as good as it is in isotropic ma- 
terials. This can be seen in a number of previous figures, for example, 
fig. 54. Thus the discrepancies between theory and experiment may result 
from not being able to accurately locate the center of a fringe. Fringe 
count was kept low in the orthotropic material, Ny - 0.5, whereas it was 
made high in the isotropic material, N * 5.0. Due to streaking by the 
fibers, the higher order fringes (greater than 3 or 4) are difficult to 
define in the orthotropic material. In an isotropic material, higher 
order fringes are narrower and more distinct than the lower order 
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Fig. 78 Comparison between predicted and observed isochromatic 
fringe order along 30° diameter, isotropic disk. 
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fringe order along 30° diameter, isotropic disk. 


152 



fringes. Thus the two materials respond in an opposite manner to high 
fringe numbers. Where the stresses are high and changing rapidly, the 
isotropic material responds in an ideal manner. The orthotropic materi- 
al responds exactly opposite to the ideal manner. This lack of spatial 
resolution of fringes in the orthotropic material is felt to be the ma- 
jor reason for the discrepancies between theory and experiment. Also, 
perhaps a lack of spatially uniform residual bi refringence contributes 
to the behavior. 

The real issue, though, is how accurately can stresses be predicted 
using orthotropic photoelasticity. This is the subject of the next 
section. 
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BACKWARD TESTING: PREDICTION OF STRESSES 


Historically, photoelasticity has been used to predict stresses 

from known isochromatic and isoclinic fringe information. For a general 

state of stress, at each point in a body the two normal stresses, 

and a , and the shear stress, t , are to be determined. At each point 
y xy 

in the body, the isochromatic fringe count gives one quantity of experi- 
mental information and the isoclinic gives a second quantity of experi- 
mental information. A third quantity or condition is needed if three 
stresses are to be uniquely determined at the point. In isotropic ma- 
terials, interferometric measurements have been used to provide informa- 
tion regarding the sum of principal stresses. This is a third quantity 

and thus a , a , and x can be determined. Interferomet ri c measure- 
x y xy 

ments are not that easy to obtain so experimenters tend to rely on other 
relations among the stresses as a third condition. One other condition 
is that the unknown stresses must satisfy the compatibility equations. 
Compatabi 1 i ty is a second order partial differential equation. For iso- 
tropic materials the compatibility equation reduces to Laplace's equa- 
tion for the quantity + a^). Since the solution to Laplace's equa- 
tion on simple domains such as a circle or rectangle is well estab- 
lished, investigators have used the compatibility equation as the auxil- 
iary condition. For orthotropic materials, the compatibility equation 
in terms of stresses does not reduce to Laplace's equation. 

Another condition often used by photoelasti cians is that the un- 
known stresses must satisfy equilibrium. Equilibrium is represented by 
two first-order partial differential equations. These two equations in- 
volve the three stresses and by themselves cannot be solved, as 
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Laplace's equation can, on any domain. In addition, the two equilibrium 
equations plus the Isochromatic and the Isoclinic fringe data result in 
four conditions from which to solve for the three stresses. In this 
situation the stresses are overdetermined. To have a unique solution, 
one of the equilibrium equations can be simply discarded. This is done 
In the so-called shear difference approach. If both equilibrium equa- 
tions are retained the stresses computed are the three stresses which 
satisfy the four conditions in a least-squares sense. It is also possi- 
ble to use the equilibrium equations in such a way that in the region of 
Interest In a photoelastic model there are exactly as many unknown 
stresses as there are equations. This approach was used in this present, 
study and it relies on applying the equilibrium equations in an average 
sense* 

Independent of how the equilibrium equations are used as an auxili- 
ary condition, the equilibrium equations are usually represented in fi- 
nite-difference form. The partial derivations are approximated by using 
central, forward, and backward differences in the stresses at the gird, 
or mesh, points of a finite-difference grid. The partial differential 
equations thus become algebraic equations relating the three stresses at 
the mesh points. In isotropic photoelasticity, eqs. 78 and 79 are then 
used as two more algebraic equations at each mesh point. Since the iso- 
chromatic, N, and the isoclinic, 0, are known at each point, eqs. 78 and 
79 represent algebraic relations the normal and shear stresses must sat- 
isfy at each point. The two finite-difference forms of the equilibrium 
equations and eqs. 78 and 79 represent the four equations which are used 
to compute the three stress components. A previous photoelastic study 
[14] of isotropic double-lapped, double-pin connectors used this four- 
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equation scheme to estimate stresses. This scheme was adapted to 
orthotropic photoelastic materials by using eqs. 86 and 87 in place of 
eqs. 78 and 79. This is possible because the form of the equilibrium 
equations, unlike the form of the compatibility equation, does not 
depend on the material properties. The finite-di fference representation 
of the plane-stress equilibrium equation remains unaltered. Appendix C 
describes the scheme. 

To determine the degree to which orthotropic photoelasticity could 
be used to predict stresses, a final series of experiments was conducted 
with the disk using the finite-difference formulation. The disk was 
loaded in diametral compression and the isochromatic and the isoclinic 
fringe data determined along various diameters of the disk. These data 
were used in the finite-difference formulation to calculate the three 
stresses along these diameters. These stresses were then compared with 
the values computed using the elasticity solution of Appendix A. Since 
the level of residual birefringence was an issue, the experimentally 
predicted stresses were computed using various values of assumed residu- 
al birefringence. The values of N R assumed were: 0, 0.07, and 0.15. 

The value of e R was assumed to be 0. For given values of Nj and 0-p, the 

experimentally computed values of a Y , a., and t would be dependent on 

x y xy 

N R . Changing N R and e R changes the right-hand side of eqs. 86 and 87. 
Physically, changes in the right hand side of the equation 'adjust' the 
fringe data so only those data due to the applied load generate numeri- 
cal values for the stresses. Obviously if Nj “ N R and = e R , the 
right hand side 'adjusts' to zero and the system of equations predicts 
that there are no stresses due to applied loads. 
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Before using actual photoelastic data in the finite-difference 

scheme, the accuracy of the finite-difference formulation was assessed. 

This was done by using the exact elasticity solution and the stress- 

optic law to generate Ny and 0^ at each point of the finite-difference 

mesh. These values of Ny and were used in the finite-di fference 

scheme as if they were experimental data. The values of cr , a , and 

x y 

x were then computed. These values of stress were compared with the 

xy 

values of stress at each mesh point as predicted by the elasticity solu- 
tion used to generate Ny and 0^. This exercise tested the accuracy of 
the finite-difference formulation. In essence, if the exact values of 
Ny and 0^ were used, could the finite-difference scheme compute the ex- 
act value of the stresses? Generally the answer to the question was a 
resounding yes. For a few cases, as will be seen, there were some dif- 
ferences between the exact values of a , a , and x from the elasticity 

x y xy J 

solution and the values from the finite-difference scheme. 

Figure 80 shows the finite-difference mesh used. In the following 
figures the stresses along the middle line of the three-line mesh, i.e. 
the diametral line, are reported. The mesh consisted of 0.002 in. by 
0.092 in. squares. With the disk 3.12 in. in diameter, there were 18 
grid points along the disk radius. Studies were not conducted to deter- 
mine the effects of mesh refinement on predicted stresses. Because of 
the aforementioned good agreement between the elasticity solution stres- 
ses and the finite-difference predictions using the exact values of Ny 

and 0y , such refinement studies were felt to be superfluous. 

Figures 81-83 show the stresses along the diameter perpendicular to 
the load for the case of the fibers aligned with the load. The stresses 
have been normalized by the absolute value of the mean stress across the 
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Fig. 80 Finite-difference mesh on disk. 
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disk diameter, |P|/Dt, where P is the load, 0 is the disk diameter, and 
t is the disk thickness. (In the elasticity solution P is negative if 
compressing the disk.) The normalized stresses are plotted on the ver- . 
tical axes and the distance along the diameter is plotted on the hori- 
zontal axis. The diametral distance is measured from the left edge of 
the disk. Experimental measurements of Ny and 0 ^. were obtained for only 
one-half the disk. The inset in each figure indicates the fiber orien- 
tation relative to the load, and the diameter being considered. 

Figure 81 shows variation of the major compressive stress, a , with 

J 

diametral distance. There are several sets of data on this and on sub- 
sequent figures. The solid squares represent the predictions of the 
elasticity solution. The open diamonds represent the predictions of the 
finite-difference scheme when the elasticity solution values of Ny and 
O j were used as 'experimental* data. This is the aforementioned finite 
difference check. In fig. 81 the diamonds are almost entirely obscured 
by the solid squares, indicating the aforementioned good accuracy of the 
finite-different approach. The slight differences that did occur were 
at the outer edge of the disk where the stresses were quite low. The 
open squares represent the stress predictions when the values of Ny and 
0-j., as measured from the actual disk, are used in the finite-difference 
scheme and the residual birefringence was assumed to be Nr = 0.15. The 
open triangles indicate the experimentally predicted stresses when Nr 
was assumed to be 0.07. The open circles indicate the case when Nr was 
assumed to be zero. As can be seen, and as could be expected, the level 
of assumed residual bi refringence affected the stress predictions. The 
close agreement between the elasticity solution (solid squares) and the 
Nr = 0.07 case was in keeping with the findings of the forward testing 
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experiments described in the previous section. There, recall, with N R = 
0.07 or lower, the level of residual birefringence resulted in better 
correlation with experimental findings than other levels of assumed 
residual bi ref ringence. 

Figure 82 shows the matrix direction stress, o x , for this same 
case. In this figure the open diamonds are not obscured by the exact 
solution's solid squares. It is important to note these stresses are an 
order of magnitude less than the stresses of the previous figure and 
the scale for the stress coordinate has been expanded to fill the fig- 
ure. Plotting this figure with the same scale as used in fig. 81 would 
significantly 'squeeze' the a data together. What was interesting 

A 

about the experimental a computation is that it did not depend on the 

A 

assumed level of residual bi ref ringence. This is evidenced by the fact 
that the open circles, squares, and triangles are coincident. This 
always occurred when the x-y system coincided with the principal materi- 
al axes. 

Figure 83 shows the shear stress calculations, t , for the case 

xy 

being discussed. Due to symmetry in geometry, loading, and material 
properties the shear stress was zero across the horizontal diameter. 

The shear stresses shown in fig. 83 were the values for the upper mesh 
points of the three-line fi nite-di fference grid of fig. 80. These shear 
stress were quite small and really almost meaningless. The calculations 
could be easily influenced by any inaccuracy in measurement or computa- 
tion. However, the trend and magnitude of the theoretically predicted 
shear stresss were well matched by the experimentally determined values. 
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Figures 84-86 show the three stress components along a diameter 
oriented 30° relative to the horizontal. The stresses are computed in a 
coordinate system which has its x axis along this diameter and the y 
direction perpendicular to this. The orientation of the x-y system is 
shown in the figures. 

Figure 84 shows the variation of ay along this diameter. The 
symbol designations are, and will be, consistent with previous notation. 
Here the stress calculations based on the experimental values of N T and 
0-j- did not seem to align with any particular level of assumed residual 
bi refringence. However, the experimental calculations were more closely 
aligned with the theoretical predictions when low levels of residual 
bi refri ngence were assumed in the finite-difference calculations. 

Figure 85 illustrates the a x stress computations. Again the scale 
differences between this figure and fig. 84 should be noted. 

Figure 86 illustrates the t stress calculations for the 30° 

xy 

diameter. These stresses were on the true diameter rather than for an 
off-set line as in fig. 83. For the majority of the diameter, the ex- 
perimental calculations using Nr = 0.07 coincided with the theoretical 
elasticity solution. 

Figures 87-104 show the stresses for other load orientations rela- 
tive to the fibers. These situations will not be discussed in detail 
except for cases where there were some unusual features. In some cases 
the correlation between the exact solution and calculations based on ex- 
perimental data was similar to the reasonable correlation for the cases 
just discussed in detail. In the other case, correlation was not as 
good. 
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Fig. 87 Comparison of a stresses along diameter perpendicular to load, 
fibers 30° relative to load. 
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Fig. 90 Comparison of a stresses along diameter perpendicular to fibers, 
fibers 30° relative to load. 
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Fig. 93 Comparison of a stresses along diameter perpendicular to load, 
fibers 45° relative to load. 
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fibers 45° relative to load. 
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Fig. 97 Comparison of a stresses along diameter perpendicular to fibers, 
fibers 45° relative to load. 
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Fig. 99 Comparison of a stresses along diameter perpendicular to load, 
fibers perpendicular to load. 
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Fig. 100 Comparison of a stresses along diameter perpendicular to load, 
fibers perpendicular to load. 
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Fig. 103 Comparison of q stresses along 30° diameter, fibers perpendic 
ular to load. 


184 



/( |P|/Dt) 


O experiment, N R assumed 0 
A experiment, N R assumed 0.07 



Fig. 104 Comparison of t stresses along 30° diameter, fibers perpendicular 
to load. xy 


185 




The shear stresses in fig. 89 showed a significant sensitivity to 
the assumed value of Nr. It appeared that if a value of Nr = 0.02 or 
0.03 had been used, the correlation between theory and experiment for 
the shear stresses for that particular load and diameter orientation 
would have been good. The a stresses for this case, fig. 87, would 

J 

have also correlated well. To keep the scatter of fig. 89 in context, 
however, the scale of fig. 89 should be compared with the scale of fig. 
87. 

Fig. 90 shows sensitivity of the major compressive stress component 
to the assumed value of Nr. Here again a value of Nr = 0.02 - 0.03 
would have given excellent correlation. 

Figure 91 is the first case to show differences between the exact 
elasticity solution and the finite-difference calculations (solid squares 
and open diamonds, respecti vely) when the exact values of Ny and were 
used. Due to the magnitude of the stress in Fig. 91 relative to the 
major stress, fig. 90, there was no concern for the differences. 

The case shown in figs. 93-95 deserves mention. Referring to eq. 
85e, the value of <J> for figs. 93-95 is <f> = 45°. This means 1/C 5 is 0, 
or C5 + ». Also, from eqs. 85a and 85b, = -C£. This is the special 

case discussed in eqs. 89a-91. It required a different computer pro- 
gram. The exact solution and the finite-difference calculations, solid 
squares and open diamonds, were in good agreement for this case. This 
gave confidence in this special -case computer program. This particular 
load/fiber orientation resulted in poor correlations between theory and 
experiment, e.g. fig. 93. 

The cases where the load was perpendicular to the fibers, figs. 99- 
104, offered some unusual comparisons. The exact elasticity solution 
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and the finite-difference calculations using the exact values for Ny 

and 0y deviated more in the a* stress than they did for any other case 

studied. Strangely enough, the o and t stresses, figs. 99, 101, 102, 

y xy 

and 104 showed good comparison between the two 'exact' cases. 

As in the previous section on forward testing, it was of interest 
to determine how closely the stresses could be computed experimentally 
in an isotropic disk. Figures 105-107 show the three components of 
stress, in an isotropic disk, across a diameter perpendicular to the 
load. The shear stress was computed for the top line of the three-line 
finite-difference grid. There was no residual birefringence in the 
material . 

Figure 105 shows the major compressive stress, o^ m Here as in the 
previous discussions, the solid squares are the exact elasticity solu- 
tion while the open circles are the results computed from the experimen- 
tal values of N and 0 and the finite-difference scheme. The comparison 
was good. Exact values of N and 0 from the elasticity solution were 
used to test the finite-difference calculations and the results of that 
calculation are displayed as open diamonds. It is impossible to see 
them in fig. 105 since they coincide with the solid squares. 

Figure 106 shows the a ^ stress calculations while fig. 107 shows 
the shear stress calculations. The comparisons between the exact solu- 
tion and the experimental calculations were good for these stresses. In 
fact, they were much better than the comparisons for the minor stresses 

(a and t ) in the orthotropic material. Figures 108-110 show the 
x xy 

stress comparisons for a diameter which was 30° relative to the diameter 
perpendicular to the load. The correlation was again good. 
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Fig, 105 Comparison of a stresses along diameter perpendicular to load, 
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From the results for the orthotropic and isotropic disks, it is 
clear, as it was in the section on forward testing, orthotropic 
photoelasticity does not appear to be as accurate as isotropic photo- 
elasticity. One is faced with the question. Is it the hypotheses and 
assumptions that are lacking, or is it the practicing of the theory? 

The lack of good fringe definition could certainly lead to scatter in 
the stresses predicted from the experimental data. Uncertainty in, or 
simple miscalculation of, the residual birefringence certainly could 
also produce the poor results. As has been seen, simply adjusting the 
level of residual birefringence can lead to excellent correlation be- 
tween the theoretical elasticity prediction of stresses and the value 
determined from the experimental photoelastic data. All evidence seemed 
to indicate the level of residual bi refri ngence, and perhaps a lack of 
spatial uniformity, were responsible for the lack of correlation that 
can be achieved with an isotropic disk. To look at the lack of correla- 
tion from one other point of view, a numerical study of the effects of 
inaccuracies in the isocline measurements was conducted. This study, 
and the motivation for it, are discussed in the next section. 
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EFFECTS OF ISOCLINE ERRORS 


ON STRESS PREDICTIONS 


From the figures and discussion of the fringe patterns in the or- 
thotropic disk, e.g. fig. 51, it is clear that fringe definition was a 
concern in the experiments. The figures of fringe patterns shown previ- 
ously are illustrations of the isochromatic fringe. The isoclinic 
fringe have not been illustrated. Several isoclinic fringes, along with 
their computer-generated predictions, could have been shown in the sec- 
tion on forward testing: global behavior. Since for stress predictions 

local behavior was more of a concern than global behavior, that exercise 
was not pursued. However, the importance of the isoclines should not be 
under-emphasised. In fact, in isotropic photoelasticity, the stress 
calculations are generally much more sensitive to inaccuracies in the 
isoclinics than to inaccuracies in the isochromatics. In isotropic pho- 
toelasticity, the isoclinic fringes are much harder to discern than the 
isochromatic fringes. The isoclinics are obscured by the isochromatic 
fringes and they are not as sharp. In addition, in regions of low 
stress gradients it is very difficult to see the isoclines. So it was 
with these experiment with orthotropic materials. The optical isocline 
was, at times, much more difficult to see than the isochromatics. Dis- 
persion by the fibers, regions over which the isoclines did not change 
rapidly, and lack of vivid contrast between light and dark regions of 
the model, all contributed to uncertainties in the isoclinic data. For 
this reason, numerical studies were conducted to determine the effect of 


uncertainties in the isocline data on the predictions of a 


t for the various fiber and diameter orientations. There were strong 

xy 
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indications that some of the discrepancies observed between the exact 
stress calculations and the computations based on experimental data were 
due to uncertainties in the residual birefringence level. However, for 
completeness, the effects of uncertainties in the isocline measurements 
were studied. 

The isocline uncertainty studies were of two types. Both were 
based on the isoclines predicted by the exact elasticity solution and 
the stress-optic laws. These were the same data used to generate the 
open diamonds in figs. 81-104. In the first uncertainty study, the the- 
oretical value of the isocline of each mesh point in the finite-differ- 
ence scheme was increased by 5°. This adjusted isocline was then used 
with the exact Ny as 'experimental data' in the finite-difference calcu- 
lations. All stresses of all eight fiber/diameter orientations were 
computed using the skewed isocline. These stresses were compared with 
the stresses from the finite-difference calculations using the correct 
values of e T and, of course, Ny. The results of this study were used to 
determine how a 5° error in all isoclines affected the predicted stres- 
ses. The value of 5° was chosen as a value which could actually occur 
if several experimental effects were additive. The polariscope being 
used was a split bench type and the polarizer and analyser were linked 
together by shafts, universals, and other coupling mechanisms. In addi- 
tion, these optical elements were rotated by means of gears. It is en- 
tirely possible that backlash in the gears and coupling elements could 
produce a constant bias in the measurements. In addition to the mechan- 
ical bias in the polariscope, the alignment of the disks in the polari- 
scope was another source for a constant bias. It is possible that the 
fiber direction of the disk was not aligned exactly as assumed in the 


196 



polari scope. There is a 0° direction in the polari scope but there were 
no cross-hairs or other alignment aids to insure that the direction of 
polarization of the optical elements were aligned at a particular angle 
relative to the disk's fibers. The polariscope was equipped with the 
standard devices for measuring isocline angle, i.e, a pointer and gradu- 
ated marks on the circumference of the analyser. The disk could have 
been rotated inplane by a degree or two relative to the zero degree dir- 
ection of the polariscope. Related to this, it was evident that within 
the disk the fibers were not perfectly straight. Thus the fiber/pol ari - 
scope alignment could have been biased in several ways. No effect was 
assumed to be large. It was simply assumed all effects added to about 
5°. 

In the second numerical study, the exact optical isocline at each 
mesh point was perturbed by adding to it an angle that was generated by 
a Gaussian random number generator. The mean of these angles added to 
the exact angles was 0° while the standard deviation was 5°. This ran- 
dom perturbation represented possible experimental error in interpreting 
the center of the isoclinic fringe. The error was assumed to result 
from the aforementioned lack of good definition of the isoclinic 
fringes. All stresses for all 8 fi ber/di ameter orientations were compu- 
ted and compared with the finite-difference computations using exact 
values of Qj and Ny of each mesh point. 

As expected, a rnyriad of results were generated from the uncertain- 
ty studies. The effects of the two types of perturbati ons varied wide- 
ly, depending on the fiber and diameter orientation, and the particular 
stress component. For example, the a stress along the diameter perpen- 

y 

dicular to the load, with the fibers aligned with the Toad, was not af- 
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fected by either type of error. The same was true of the shear stress 
for the case with the fibers 30° relative to the load and the diameter 
perpendicular to the fibers. On the other hand, there were cases where 
one type of error had an effect but the other type did not. For exam- 
ple, the stresses along a diameter perpendicular to the load and with 
the load perpendicular to the fibers were more sensitive to random er- 
rors than to the constant error. Finally, some stresses were seriously 
affected by either disturbance to the isocline. The shear stresses 
along a diameter perpendicular to the load and with the fibers 30° rela- 
tive to the load were very sensitive to isocline errors. 

Tables 4 and 5 rank the various situations as to how they were af- 
fected by the errors. The top of the column corresponds to conditions 
least affected. There is no significance to the fact that a is the 

y 

left column, a the middle column and t is the right column. That is, 
x xy 

it cannot be said that the least affected a stress was less affected 

y 

than the least affected a* stress. It is possible to say, however, that 
for a specific fiber orientation and diameter, the larger of the three 
stresses was generally least sensitive to errors. In the tables the 
figure number which illustrates the fiber orietation relative to the 
load, and the particular diameter, is indicated in parenthesis. Obvi- 
ously these figures do not represent the perturbed stress computations. 

The major issue, however, was: Could the discrepancies in the 

stresses between theory and experiment observed in the various situa- 
tions for the actual disk be explained by the presence of either a con- 
stant or random error? Figure 111 shows the predicted stresses along 
the diameter perpendicular to the fibers for the case when the fibers 
were 30° relative to the load. The figure shows the values predicted 
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Table 4 


Ranking of Effects of Constant Isocline 
Error on Stress Predictions ^ 


a 

y 


fiber/load diameter 


angle 

(fig. no.) 

ali gned 

perp. to 
load(81) 

perp. 

perp. to 
load(99) 

aligned 

30° (84) 

perp. 

30° (102) 

45° 

perp. to 
1 oad(93) 

30° 

perp. to 
load(87) 

30° 

perp. to 
fibers (90) 

45° 

perp. to 
fibers (96) 


a 

x 


fi ber/load 
angle 

di araeter 

al i gned 

perp. to 
load(82) 

perp. 

perp. to 
load(100) 

45° 

perp. to 
f1bers(97) 

45 

perp. to 
load(94) 

aligned 

30° (85) 

30° 

perp. to 
fibers(91) 

30° 

perp. to 
load(88) 

perp. 

30° (103) 


T 

xy 

fiber/load diameter 
angle 


30° 

perp to 
fibers (92) 

al igned 

30° (85) 

45° 

perp. to 
fibers (98) 

al igned 

perp. to 
load(83) 

perp. 

30“ (104) 

perp. 

perp. to 
load(102) 

45° 

perp. to 
load(95) 

30° 

perp. to 
load(89) 


(!) least affected at top 
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Table 5 


Ranking of Effects of Random Isocline 
Error on Stress Predictions^*) 


0 0 T 

y x xy 


fiber /load 
angle 

diameter 
(fig. no.) 

fi ber/1 oad 
angle 

diameter 

fi ber/1 oad 
angle 

diameter 

a 1 i gned 

perp. to 
load(81) 

45° 

perp. to 
fi bers (97) 

30° 

perp. to 
fi bers(92) 

perp. 

perp. to 
load(99) 

30° 

perp. to 
fibers(91) 

al i gned 

30° (85) 

aligned 

30° (84) 

a 1 i gned 

30° (85) 

45° 

perp. to 
fibers (98) 

0 

If) 

perp. to 
1 oad(93) 

30° 

perp. to 
load(88) 

aligned 

perp. to 
1 oa d( 83 ) 

perp. 

30° (102) 

perp. 

30° (103) 

perp. 

30° (104) 

30° 

perp. to 
load(87) 

al i gned 

perp. to 
load(82) 

perp. 

perp. to 
load(102) 

30° 

perp. to 
fibers (90) 

perp. 

perp. to 
load(100) 

30° 

perp. to 
load(89) 

45° 

perp. to 
fibers(96) 

45° 

perp. to 
1 oa d (94) 

45° 

perp. to 
load(95) 


(*) least affected at top 
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using the exact values of the isocline and predictions based on the two 
types of disturbances to the isoclines. This figure should be compared 
to fig. 91. The case of the constant error in fig. Ill looks somewhat 
like the mismatches in fig. 91. Figure 112 shows the predicted o 
stresses along a diameter perpendicular to the fibers for the case of 
the fibers making a 45° degree angle relative to the load. This figure 
should be compared with fig. 96. 

From the figures it is clear a constant error was much more likely 
than random errors to explain the discrepancies observed between theo- 
retical and experimentally measured stresses. However, there was a sig- 
nificant difference between the discrepancies observed in the actual 
disk and the discrepancies generated in the numerical error studies. In 
both fig. Ill and 112 it is seen that the constant isocline errors 
seemed to affect stresses near the center of the disk more than they af- 
fected the stresses near the disk's edge. For example, in fig. 112, the 
triangles and circles coincide up to about x/D = 0.15. In fig. Ill 
these symbols coincide up to x/D = 0.2. After those locations, the two 
predictions diverge as the center of the disk is approached. In con- 
trast, fig. 96, for example, shows a constant discrepancy along the 
whole diameter, starting at the edge of the disk. This would seem to 
indicate a constant error of a type different than a constant isocline 
error. As was felt all along, and as these isocline error studies con- 
firmed, errors in the level of residual bi refri ngence were more likely 
to account for discrepancies than other types of errors. 

For completeness, similar scatter studies were conducted for the 
isotropic disk stresses. Recall, the experimentally determined stresses 
for the isotropic disk were quite close to the predictions of the 
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Fig. 112 Effect of isocline errors on predicted a stresses along 

diameter perpendicular to fibers, fibers y 45° relative to load. 
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elasticity solutions. Figure 109 illustrates the stress which deviated 
the most from predictions. Figure 113 shows the stresses computed for 
that same situation but with the two types of errors incorporated into 
the isocline. As can be seen, both types of errors affect this stress 
calculation to a significant degree. Comparing figs. 109 and 113, it 
appears something other than random errors caused deviation in fig. 109. 
However since neither type of error caused the constant shift along the 
entire diameter, as for example in fig. 89, the isotropic disk was not 
studied further. 

The results of the error study indicated that while isocline errors 
can lead to serious errors in the prediction of stresses, such errors 
were not likely the source of the errors observed in the experiments 
with the orthotropic disk. The errors observed had very few of the 
characteristics of random errors and they did not seem to behave exactly 
like a constant error. These errors were thus discounted as a major 
source of discrepancy. 
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SUMMARY AND CONCLUDING COMMENTS 


A considerable amount of information, both in the form of the deri- 
vation of equations and in the form of data, has been presented. In the 
second section equations governing the photoelastic phenomenon in ortho- 
tropic materials were derived. In the third section, calibration of the 
material was discussed and calibration data were presented. The fourth 
and fifth sections dealt with testing and checking the orthotropic pho- 
toelastic theory in the context of a unidirectional disk in diametral 
compression. The just-presented section discussed sources of errors in 
photoelastic materials. As data were presented, they were discussed and 
put into context with other data. The purpose of this section is to 
summarize the important results presented. 

The theory of orthotropic photoelasticity was developed in the sec- 
ond section so that: (1) the effects of residual birefringence could be 

accounted for in the photoelastic phenomenon, and (2), the key steps 
leading to specific operational equations, e.g. stress prediction equa- 
tions, were shown. The experiments for calibrating the photoelastic ma- 
terial revealed that uniaxial tests are not the best way to calibrate 
the material. The four-point bend specimen was found to be superior to 
the uniaxial specimen for calibration purposes. The reasons were as 
follows: 

(1) The stress gradients in the bending specimen will lead to 
sharper fringes and thus better fringe definition. 

(2) A single fixture, specimen, and test can be used to load the 
material in both tension and compression. 
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(3) The optical neutral axis is a convenient way to assess residual 
birefringence. 

m- 

Unfortunately, uniaxial tests must be conducted for determining elastic 
properties and so the four-point bend tests constitute an extra set of 
tests. 

The comparisons of the computer-generated and the observed fringe 
patterns revealed that on-a-whole, the notion of orthotropic photoelas- 
ticity, as presented, is valid. As with the uniaxial tests, sharpness 
of fringes was shown to be a concern. Also, it was shown that the level 
of residual birefringence affected the correlation between the predicted 
and observed fringe patterns. By adjusting the assumed level of residu- 
al birefringence the predicted and observed fringe pattern could be made 
to coincide for a particular case. Unfortunately, the adjustment was a 
function of the load orientation relative to the fibers. The computer- 
generated fringes also showed that in the presence of residual birefrin- 
gence, the sign of the applied load can strongly affect the fringe pat- 
tern. 

A detailed examination of the predicted and observed photoelastic 
behavior was provided by studying the isochromatic and isoclinic fringes 
along specific disk diameters. This type of point-by-point analysis was 
necessary for accurate determination of the stresses. The fringe pre- 
diction study was contrasted by conducting similar studies on isotropic 
disks. From the fringe prediction studies, it can be concluded that: 
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1) The material seemed to exhibit behavior associated with lower 
rather than higher levels of residual birefringence. Specific- 
ally, a value between 0 and 0.07 for N R would have correlated 
well with 80% of the cases. 

2) Compared to isotropic photoelasticity, orthotropic photoelas- 

ticity, as hypothesised here, does not seem as accurate. This 
may be related to the lack of good fringe resolution or more 
basically, to the assumptions of the theory itself. Three as- 
sumptions which could be suspect are: (a) validity of through- 

the-thickness averaging, (b) validity of the stress-optic law 
being orthotropic, and (c), validity of the superposition of 
residual and mechanically-induced birefringent effects. To 
challenge these assumptions requires basic research beyond the 
scope of this study. 

The section on stress prediction reinforced the idea that residual 
birefringence in the material was less than N R = 0.15. Specifically a 
value between 0.0 and 0.07 appeared to be correct. It is very important 
to note, however, in no computation nor experimental measurement was 
there overwhelming support for the assumption N R = 0. Some value of N R 
could be chosen so that correlation was good. Unfortunately, the level 
of N r needed for good correlation varied from one situation to the 
other. However, it can be concluded that, in general, stresses can be 
predicted with confidence using orthotropic photoelastic material. Un- 
fortunately in stress prediction, and in every other phase of the corre- 
lation between theory and experiment, when the fibers were perpendicular 
to the load the correlation between theory and experiment was poor. 
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Finally, studies examining the effects of isocline errors showed 
that while isocline errors can lead to serious consequences, they were 
not likely to have caused the observed discrepancies between the 
experiments and the theoretical predictions. 
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APPENDIX A 


The form of the dielectric tensor 


The dielectric tensor can be expressed as 


/ e ' 

1 / 1 + x \ 

\ X 1 

1 \ * 

) e < 
y 1 

II 

H 1 
+ 

5> 

<U 



where x v » Xw> and X yv are components of the susceptibility tensor . 

* y 

Numerically, the components of the susceptibility tensor are small 
compared to unity. The quantity e aye is defined as 


e + e 
J< X 


ave 


(A— 2 ) 


Substituting eq. A-l it is seen that 


e ave 


1 + 


x x + x y 

2 


and thus e ave is of order unity. 
The quantity Ae is defined as 


* -V * *5, • 

Using eq. A-l, Ae becomes 


(A-3) 


(A— 4) 
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AS = z y ~) ♦ x^y • (A-5) 

The quantity Ae is small compared to unity and in particular, it is 
small compared to e ave . Thus the quantity 



+ Ae 


(A-6) 


can be approximated, using the binomial expansion, as 


Similarly, 


(A— 7) 


e ave 



1 Ae 

2 e 


•) • 


ave 


( A -8 ) 
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APPENDIX B 


Elasticity Solution for Stresses in an Orthotropic Disk 
in Diametral Compression 


The stresses in an orthotropic disk in diametral compression can be 
represented by a state of plane stress. If the x-axis is aligned with 
the stiff (fiber) direction and the y-axis is aligned with the soft (ma- 
trix) direction, then the z-axis is perpendicular to the plane of the 
disk. The nonzero stresses in the disk are a x , o^ 9 and t . (In the 
discussion here, x is used to denote fiber direction. In the text this 
direction has been denoted as 1M. Likewise the y direction here is 
denoted 2M in the text. Using x and y instead of 1M and 2M is done 
simply so the subscripts on a, e, y, x, etc. are more familiar.) 

The three-dimension equilibrium equations, in rectangular cartesian 


coordinates, reduce to 


and 


da x d V 

— - + — = 0 , 

ax ay 


dT xv aa y 

-IL + = o. 

ax ay 


Hooke's Law becomes 




(B-l ) 
(B-2) 

(8-3) 


The only nontrivial compatibility equation is 


a 2 y 

a 2 e . 

a 2 e.. 

.. JJ& = 

* + -i 

ax ay 

ay 2 

ax 2 


( 8 - 4 ) 
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Typically in plane stress problems, a stress function F(x,y) is defined 
such that 


_ 

x ay 2 

(B-5a) 

a 2 F 
" ax 2 

(B-5b) 

a 2 F 

T xy “ “ axay • 

(B-5c) 

With eq. B-5, equilibrium is automatically satisfied. 

Using Hooke's 


Law, eq. B-3, and eq. B-5, the compatibility equation becomes: 


a™ + ( 2 a, „ + a 


22 


ax' 


12 


66 
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d 4 F 


ax 2 ay 2 


- + 


a,, ^ F - - 0 

11 ay 4 


( 8 - 6 ) 


Equation B -6 admits a solution of the form 

F = F(x + py) . (8-7) 

where p is a constant. 

Substituting eq. B-7 into eq. B -6 results in the so-called 
characteri Stic equation, 

a ll^ 4 + C 2a 1 1 + a 66^ 2 + a 22 " 0 • (8-8) 


The solution given by eq. B-7 is only a solution if p assumes spe- 
cific values, namely one of the 4 roots of eq. B- 8 . The roots of eq. 
B -8 are generally complex, the roots being denoted p^, p 2 * P 3 and p^. 
The complete solution of eq. B -6 is thus 


F = F x (x + (ijy) + F 2 (x + (J 2 y) + F 3 (x + n 3 y) + F 4 (x + ^y) . (B-9) 
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The roots are in pairs, two roots being complex conjugates of the other 
two, i.e. 

P 3 = Pj and p 4 = p 2 • (B-10 a,b) 

It is convenient, and the basis of the complex variable approach, 
to define 4 complex variables z^, z 2 , Z 3 , and z 4 as follows: 


z 1 = x + , z 2 = x + , z 3 = x + ^y and z 4 = x + p 4 y . 

(5-11 a,b,c,d) 

Obviously by eq. B-10, 

z x = and - z^ . (B-12 a,b) 

By eqs. B-ll and B-12, eq. B-9 can be written as 

F = F 1 (z 1 ) + F 2 (z 2 ) + + F 4 (i 2 ) . (B-13) 

Since the stresses are real, F must be real and so eq. B-13 can be 
written as 


F = aRefF^Zj) + F 2 (z 2 )] . 


(B-14) 


To use eq. B-14 in eq. B-5, the following partial derivatives of F 
are needed: 

ic dF, dF, 

(8-15a) 


dL = d A + d !i 

5 x dZ| dz 2 


eL . d . F _i + 
dy ^1 d~z 1 **2 dz 2 


(3-1 5b ) 
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II 


Usually the derivatives on the right hand sides of eq. B-15 are 


redefined, i.e. 

d F,(z.) 

♦l (z P - 

{ 8- 16a ) 


dF (z 2 ) 



<t> 2 ( z 2 ) = ~ i iz~ • 

(B-16b) 

then 

— F - = *j(z.) + *i(z ? ) 
ax 2 1 1 ‘ d 

(B-17a) 


AW + 

(B-17b) 


b 2 F 

ay ax' = H ^ z l ^ + ’ 

(B— 17c ) 

where 

d<t>, (z. ) 

= - 4 -— 

(B-18a) 

and 

d^(z 2 ) 

^ z Z> - dz 2 

(B-18b ) 

With this notation 




= 2Re[i^4>{(z 1 ) + ] * 

(B-19a) 

°y 

= 2Re[(j)^(z2) + 4>^> ( 2 2 ) -J * 

(B-19b) 

and 



T xy 

= - 2Re[p^(zj) + • 

(B-19c) 


The complex variable approach to these plane stress problems is 
convenient because (1), since it is assumed F is an analytic function at 
all but a few points in the domain, the functional form of F is known, 
and (2), the boundary conditions can be related to F through rather sim- 
ple integrals. Concerning the latter, let X n and Y n represent the x- 
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and y-direction surface tractions, respecti vely, on the boundary of the 
body, in this case the disk. 

Assume that on the disk the normal surface traction, a r , and the 
tangential surface traction, T p0 , are specified functions of 0. This 
idea is shown in fig. B-l. The quantities a p (e) and T r0 (9) cannot be 
completely arbitrary in that the disk must be in translational and 
rotational equilibrium as a rigid body. With X n denoting the x- 
component of the surface traction and Y n denoting the y-component of 
surface traction, then 

x n = x n (e) = (a r (e) cose - T re (e) sin{e)) (B-20a) 


Y n = Y n (e) = (a r (e) sine + -r re (e) cose) . (B-20b) 

The boundary conditions on F become 

0 

+ <t> 2 (z 2 ) ] = - J Y n (e)Rde + Cj (B-21a) 
e 

2Re [ ( z x ) + n 2 $ 2 (z 2 ) ] = + / X n (e)Rde + c 2 , (B-21b) 


i 0 

where and z 2 are evaluated at z = Re , i.e. the boundary of the 
disk. 

Lekhnitskii* states the solutions to ^(zj) and 4> 2 ( z 2 ) f° r this 
geometry are: 

co 

4>i ( z i ) = + ^l z l + ^ ^m^l (3-22a) 

m=2 m 


d) 2 (z 2 ) - B 0 + B t z 2 + l B m P 2 (z 2 ) . ( B— 22b ) 

m=2 m 


♦Lekhnitskii, S. G., Anisotropic Plates , translated from 2nd ed. by S. 

W. Tsai and T. Cheron, Gordon and Breach, Science Publishers, New York, 
1968, eq. 33.6. 
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The Pin,^) and P 2m < z 2 > are given by 


P lm< z l> 





R 2 (l+H 2 ) ) m } 


(B- 2 3a) 




-1 

R m ci-i^r 



R 2 (l + m|) ) m + 



R 2 (l+n 2 )) m }. 


(B-23b) 


That these are solutions to ^(z^) ancJ ^ z ?) 1S not °b v i° us to even the 
serious observer of mathematics, let alone the casual observer. These 
solutions to 4>^ ( ) and ^(z?) , and other solutions for other domain 
geometries, is what has made a mark for Letkhnitskii and others working 
in this field. 

Since ^(z^) and ^(zj?) are expressed in terms of a series, to 
match the boundary conditions given by eq. 8-21, X n and Y n of eq. B-20 
need to be expressed as a series, specifically a Fourier series in e. 
This step is universal in the complex variable approach to elasticity, 
i.e. the expanding of the known boundary data into a Fourier series. 

Assuming a p (0) and T re (0) are known, the boundary integrals on the 
right hand side of eq. B-21 can be represented by 

- / V (e)Rde + c, = a + l (a a m + a a' m ) (B- 2 4a) 

o " 1 0 m=l m m 

/ 9 x n (e)Rde + c 2 = p 0 + l (i v* + V") , (B-24b) 

o m=l 

where 
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a = e 


(B-25) 


i 0 


The quantities c^, o^, p m , are known coefficients which can be 
determined by typical Fourier analysis once the functional forms of 
X n (e) and Y n (e) are specified. Substituting eq. B-24 into the right 
hand side of eq. 8-21, and eqs. 8-22 and B-23 into the left-hand side of 
eq. B-21, and equating coefficients of like powers of o results in 4 
coupled linear algebraic equations for the A m and B m of eq. B-22. These 
equations are: 

A m + B m + of + B = - <l 
m m ml m c m 

' i l A m + ^ B m + h l l + h. l 2 B m = " Prti 
+ + 

1*1*1*111 + ti 2 t 2 B m + Mn + ^ S m = * S rn 
m - 2, 3, 4. . . . . 

The quantities t^ and t 2 are given by 

1 + in. 1 + in ? 

*1 = 1 - and *2 = 1-1^ * (B-27) 

For the case m = 1, matching coefficients in leads to 

®1 + «i 

+ ^1 + ^1 + ^1 " R (B-28a) 
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o A 

+ ^2 B l + \\ * ^2 = R 

(8-28b) 

^ A 1 + "2 B 1 + ) 2B 1 = " Ri ’ 

(B-28c) 

Equations B28a-c represent 3 equations for 4 unknowns, 
and Bj. An indeterminant solution is not a problem if the 
examined in detail. Consider only the linear terms in the 
for ^(Zj) and 4>2 ( z 2 ^ ’ n e< l* Substituting only these 
into the expressions for stress, eq. B-19 results in 

Ait Aj, 8^ f 

situation is 
expressions 
linear terms 

a x = ZRe^Aj + 

(B-29a) 

o y = 2Re[Aj + Bj 

(B-29b) 

x xy = -2Re[|ijAj + |i 2 ®i ] • 

(B-29c) 

Eliminating the taking of the real part, these three expressions can be 
written as 

a x = + ^2 B 1 + (^l) 2A l ¥ (^2 ) 2 ^1 1 (B-30a ) 

°y = t A l + \ * Bi + 8 X ] 

(B-30b) 

T xy = " [^1 A 1 + + ^1 A 1 + »*2®i ] • 

(8-30c) 


The right hand side of eqs. B-30a,b,c are related directly to eqs. B- 
28c, a, b. Thus, the leading terms in the expressions for stress, eq. B- 
30a, b,c, are 
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(B-31a) 


\ 

l 


°x = 


Pi ' Pi 

Ri ” 


“1 + “1 
a v ~ R 


Pi + Pi 

T xy R 


(B-31b) 
( B— 31c ) 


Since only the derivatives of ^(zj) and <i>2 ( z 2 ) are rec I u1 * recl for the 

computation of stress, A Q and B Q in eq. B-22 do not enter into the calc- 

ulations. Thus eqs. B-31a-c and the solution to Eq. B-26 for A m and B m> 
m > 2, give the necessary constants to completely determine the state of 
stress in the disk. Evaluating ^(z^ and ^(z^) at a specific x and y 
(or r and 0 ) gives the stresses at a specific location in the disk. 

For the diametral concentrated load, the constants 0 ^ and p m , 
m > 0, are easily determined. As shown in fig. B-2, consider the 
concentrated load P at an angle a relative to the x-axis. For this 
loading 

x n (e)Rde = [P6(e - a) + P6(e - a - n) ] COSGRde (B-32a ) 

Y n (e)Rde = [P6(0 - a) + P6(e - a - n) ] sinGRdG, (B-32b) 


0 < 0 < 2n 


where 6( ) is the 'Dirac delta and +P represents diametral tension. The 
integrals on the right hand side of eqs. B-21 are 1/2-cycles of square 
waves, i.e. 
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I 


/ X n ( 0 )Rd 0 - P[U (0 - a) COSa + U ( 0 - a - it) CO$(a - tc) ]R 
0 

(B-33a) 


e 

/ Y n (0)Rd0 = P [ U ( 0 - a) siria + U(0 - a - it) sin(a - ix) ]R , 
o 

(B-33b) 

U( ) denoting the unit step function. Since 


COS ( a - n) = -COSa and - -sina, 

(B-34a,b) 

eq. B-33 becomes 


0 

/ X n (0)Rd0 = Pcosa[U{ 0-a) - U( 0 - a - tc) ]R (B-35a) 
0 


0 

/ Y n (0)Rd0 = Psi na[li( 0 - a) - U( 0 - a - n) ]R . (B-35b) 
o 


These integrals are shown graphically in fig. B-3. Using the usual 
ideas associated with Fourier analysis, the two integrals can easily be 


expanded to a 

Fourier series of the form eq. B-24. The 

results are: 


Psi na . Q Pcosa 

a o = - 2 “ and p o ~r~ • 

(B-36a,b) 

For m odd: 

a m = ( sin ( ma ) + 1 cos(ma)} 

(3-36c) 


P m = - {sin(ma) + i cos(ma)} 

(B-36d) 

for m even: 

II 

O 

II 

• 

(3-36e,f ) 
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Obviously only a finite number of terms are actually used in eqs. 
B-22 and B-24. Concentrated loads pose the most difficult convergence 
problems for elasticity solutions, particularly if infinite series are 
truncated. In reality, any finite number of terms results in a conver- 
gence problem right at the point of application of the load. So the 
question is, at points of interest removed from the point of application 
of the load, what are the convergence characteri sties? Numerical stud- 
ies, based on a FORTRAN progam which coded eqs. B-19, B-22, B-23, B-26, 
B-27, B-31, and B-36, indicated that toward the center of the disk, 10 
terms in the series gave answers nearly identical to the series with 20 
terms. Twenty terms gave answers nearly identical to a 40-term series 
over 85% of the disk. A greater number of terms caused underflow in the 
computer being used. 

As a check, a finite element program, using 216 constant-strain el- 
ements in a complete disk, was used to compute the stresses. The ele- 
ment arrangement is shown in Fig. B-4. There was no particular reason 
for using 216 elements other than the ability to be able to compute 
stresses at a large number of points. The particular finite element 
program being used computed stresses at the centroids of the elements 
and thus many elements were used. To compare the stresses generated by 
the complex variable solution with the stresses generated by the finite 
element solution, the following scheme as adopted. For a given load 
direction relative to the fibers, a, the stresses were compared along 
the loaded diameter, along a diameter perpendicular to the loaded diame- 
ter, along a diameter in the fiber direction, and along a diameter in 
the matrix direction. Such a comparison would indicate how well the 
stresses were being computed in a highly loaded region, i.e. moving from 
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the center of the disk toward the point load along a diameter, how well 
the stresses were computed in a lightly loaded region, i.e. moving from 
the center of the disk toward a free edge along a diameter perpendicular 
to the load, and whether the stress diffusion in the stiff and soft dir- 
ections was properly represented. The elastic properties of the cases 
studied were as follows: 

E x = 28.8 GPa (4.18 x 10 6 psi) 

E 2 = 9.40 GPa (1.36 x 10 6 psi) 

G 12 = 3.17 GPa (0.460 x 10 6 psi) 
v^2 = 0-3 . 

These values led to values for p (refer to eq. B-10) of: 

p x = 2.83 i , p 2 = 0.618 i . (3-37a,b) 

In the two computational schemes the disk diameter, 0, was 3.00 in. 
while its thickness, t, was unity. The comparisons are presented in 
nondimensional form, the stresses being divided by P/Dt, where P is the 
applied diametral load. The stresses are presented as a function of 
their x and y locations. The x and y coordinates are nondimensional i zed 
by dividing them by the disk radius. With the particular finite ele- 
ment, the particular mesh, and the way a point load was applied to the 
model, stresses were not computed exactly on the loaded diameter, nor on 
the diameter perpendicular to the loaded diameter, etc. For both the 
finite element and complex variable solutions, the stresses were compu- 
ted at the centroids of the elements and the element centroids were not 
exactly on the loaded diameter. 
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Tables 8-1 through B-4 indicate the comparison for 4 orientations 
of the load relative to the fiber direction, i.e. a = 0°, a = 30°, 
a = 60°, and a = 90°. Twenty terms were used in the complex variable 
solution. 

From the tables it is obvious there is good agreement between the 
finite element solution and the complex variable solution over much of 
the disk. Which is the ‘correct' solution is not clear. The finite el- 
ement mesh should be finer near the point load to represent the singular 
behavior of the stresses in that region. The singular solution is 
'built-in' to the complex variable solution but it is not clear 20 terms 
represents this portion of the solution. The problem is not important 
because in the actual experiment the disk material crazed badly near the 
point where the load was applied. Fringes could not be observed and 
compared with theoretical predictions. In addition, the load obviously 
was introduced over a finite region rather than at a point, as the 
theory assumes. Thus any fringes in the immediate region of the load 
would be strongly affected by this lack of idealism. Any comparison be- 
tween theory and experiment in this region of the disk would be meaning- 
less. 

Based on the comparison shown in Tables B-l through B-4, the com- 
plex variable solution was felt to be quite accurate and, compared to 
finite element solutions, quite quick and inexpensive to run on the com- 
puter. 
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perp. to load, 

matrix dir. load and fioer dir 


Table B-l 


Comparison of Finite Element (FE) and Complex Variable (CV) 
Solutions for Disk with Fibers Parallel to Load 
(a = 0, see fig. B-4) 


o 

c 

4-> 

c: 

Q) 




°x 

P/Dt 



P-X 


jcy 

P/Dt 


07 

X 

r 










R 

R 

FE 


CV 

FE 


CV 

FE 

CV 

i 

72 

0.100 

0.220 

0.007 

0.020 

3.26 

3.38 


3.27 

3.42 

-0.354 

-0.374 


-0.358 

-0.364 

-0.648xl0 3 

0.282x10“* 

-0.575x10“? 

0.409xl0' J 

108 

0.420 

0.040 

3.91 


4.00 

-0.376 


-0.375 

0.236 

0.202 

144 

0.650 

0.053 

5.12 


5.38 

-0.276 


-0.337 

0.788 

0.729 

180 

0.820 

0.073 

6.39 


5.46 

0.191 


0.259 

1.60 

2.14 

216 

0.950 

0.080 

7.28 


1.22 

1.72 


2.04 

2.14 

-3.34 


28 

0.007 

0.100 

2.79 

2.93 

-0.369 

-0.317 

-0.384xl0" 3 

-0.486xl0“ 3 

64 

0.020 

0.220 

2.05 

2.09 

-0.197 

-0.208 

-0.021 

-0.021 

100 

0.040 

0.420 

0.882 

0.842 

-0.048 , 

-0.055 , 

-0.017 , 

-0.021 , 

136 

0.060 

0.650 

0.198 . 
-0.279x10'} 
-0.822x10' 1 

0.168 . 
-0.250x10"} 
-0.495X10" 1 

0.476X10"- 3 

0.879X10" 13 

-0.785x10", 

-0.537x10* 

172 

0.073 

0.820 

0.011 ^ 

0.010 

0.249x10"; 

0.300xl0- 2 

208 

0.080 

0.950 

0. 369x10' J 

-0.059 

0.654x10"* 

0.017 



Table B-2 


Comparison of Finite Element (FE) and Complex Variable (CV) 
Solutions for Disk with Fibers 30° to Load 



• 





(a = 30. 

see fig. B-4) 






o 

c 

c 




a x 

P7Dt 


a 

JL 

P/Dt 



pTBt 



I 

X 

£ 










'a) 

R 

R 

FE 


CV 

FE 

CV 

FE 


CV 


33 

0.080 

0.060 

0.850 


0.887 

0.169 

0.159 

1.17 


1.17 

c_ 

69 

0.187 

0.127 

0.940 


1.00 

0.167 

0.179 

1.20 


1.20 

■5 

105 

0.347 

0.240 

1.39 


1.48 

0.209 

0.236 

1.31. 


1.33 

-a 

141 

0.533 

0.373 

3.40 


3.28 

0.251 

0.313 

1.46 


1.72 

<o 

o 

177 

0.673 

0.473 

10.6 


9.08 

0.259 

0.245 

1.29 


1.56 


213 

0.773 

0.540 

10.5 


8.58 

0.63 

0.328 

1.04 


2.08 

. 

36 

0.100 

0.007 

0.856 


0.844 

0.170 

0.169 

1.16 


1.16 


72 

0.220 

0.020 

0.822 


0.802 

0.220 

0.220 

1.17 


1.17 

■5 

108 

0.420 

0.040 

0.637 


0.643 

0.461 

0.422 

1.18 


1.17 

t- 

4) 

144 

0.647 

0.053 

0.267 


0.259 

0.917 

0.891 

1.01 


0.991 

_<D 

180 

0.820 

0.073 

-0.050 


-0.095 

1.06 

1.15 

0.445 


0.433 

«*- 

216 

0.947 

0.080 

0.016 


-0.489 

0.672 

0.684 

0.321x10 J 

-0.220 

L. 

<T3 

24 

-0.060 

0.080 

0.803 


0.878 

0.132 

0.152 

1.15 


1.14 

"a -o 

o io 

60 

-0.127 

0.187 

0.929 


0.969 

0.113 

0.129 

1.02 


1.03 

•- O 

TJ t — 

96 

-0.240 

0.347 

0.986 


1.02 

0.105 

0.086 

0.669 


0.719 

C 

41 O 

132 

-0.373 

0.533 

0.669 


0.623 

0.062 

0.054 

0.308 


0.280 

CL-tJ 

t_ 

168 

-0.473 

0.673 

0.216 


0.175 

0.014 

0.013 

0.066 


0.048 

0) 

Cl 

204 

-0.540 

0.773 

-0.074 


0.041 

-0.038 

-0.131 

-0.057 


-0.026 

C 

27 

-0.007 

0.100 

0.935 


0.921 

0.160 

0.142 

1.17 


1.14 

o 

63 

-0.020 

0.220 

1.22 


1.25 

0.073 , 

0.081 

1.04 


1.02 

X +-> 

•1- u 

99 

-0.040 

0.420 

1.38 


1.54 

-0.359x10’^ 

-0.442xl0 _Z 

-0.020 

0.614 


0.636 

«- 

4-> V 

135 

-0.053 

0.647 

0.857 


0.781 

-0.302xl0~ Z 

0.101 


0.087 

<T5 *r- 
F T3 

171 

-0.073 

0.820 

0.101 


0.058 

0.018 . 

0.025 

-0.018 


-0.014 


207 

-0.080 

0.667 

-0.198 


-0.220 

0.892xl0" z 

0.056 

-0.023 


-0.063 



matrix perpendicular fiber 

direction to load direction load direction 


Table B-3 


Comparison of Finite Element (FE) and Complex Variable (CV) 
Solutions for Disk with Fibers 60° to Load 


o 





(a = 60 , 

see fig. 

B - 4 ) 



c 




a 



0 

T 


4 -> 

c 

<v 

Er 




P70t 



P7fo 



*\) 

X 

J 

t 

FE 


CV 

FE 

CV 

FE 

CV 

30 

0.040 

0.093 

- 0.161 


- 0.157 

1.02 

1.02 

1.02 

1.02 

66 

0.093 

0.200 

- 0.166 


- 0.133 

1.06 

1.08 

1.15 

1.06 

102 

0.180 

0.387 

- 0.026 


- 0.005 

1.25 

1.28 

1.22 

1.22 

138 

0.273 

0.587 

0.535 


0.515 

1.69 

1.73 

1.73 

1.72 

174 

0.357 

0.747 

3.04 


2.82 

2.12 

2.09 

2.82 

2.95 

210 

0.400 

0.860 

8.68 


12.43 

2.74 

1.09 

3.40 

2.12 

36 

0.100 

0.007 

- 0.165 


- 0.173 

1.04 

1.03 

1.01 

0.991 

72 

0.220 

0.020 

- 0.212 


- 0.204 

1.11 

1.08 

0.913 

0.908 

108 

0.420 

0.040 

- 0.248 


- 0.254 

1.13 

1.16 

0.613 

0.623 

144 

0.647 

0.053 

- 0.180 


- 0.190 

0.913 

0.926 

0.213 

0.225 

180 

0.820 

0.073 

- 0.080 


- 0.046 

0.601 

0.597 

0.030 

0.047 

216 

0.947 

0.080 

- 0.015 


0.278 

0.381 

0.388 

- 0.022 

0.058 

21 

- 0.093 

0.040 

- 0.188 


- 0.162 

0.989 

1.00 

1.00 

0.989 

57 

- 0.200 

0.093 

- 0.162 


- 0.147 

0.964 

0.970 

0.905 

0.908 

93 

- 0.387 

0.180 

- 0.059 


- 0.089 

0.814 

0.826 

0.662 

0.674 

129 

- 0.587 

0.273 

0.002 


- 0.002 

0.574 

0.564 

0.401 

0.402 

165 

- 0.747 

0.347 

0.043 


0.074 

0.352 

0.342 

0.223 

0.224 

201 

- 0.860 

0.400 

0.053 


- 0.103 

0.218 

0.168 

0.119 

0.191 

27 

- 0.007 

0.100 

- 0.155 


- 0.146 

1.02 

1.00 

1.02 

1.02 

63 

- 0.020 

0.220 

- 0.063 


- 0.069 

0.974 

0.952 

1.05 

1.04 

99 

- 0.040 

0.420 

0.323 


0.246 

0.778 

0.768 

1.11 

1.10 

135 

- 0.053 

0.647 

1.04 


1.08 

0.361 

0.375 

0.956 

0.990 

171 

- 0.073 

0.82 

1.12 


1.31 

0.066 

0.041 

0.321 

0.328 

207 

- 0.08 

0.947 

0.073 


0.051 

- 0.019 

0.003 

- 0.052 

- 0.266 



Table B-4 


Comparison of Finite Element (FE) and Complex Variable (CV) 
Solutions for Disk with Fibers 90° to Load 
(a = 90°, see fig. B-4) 
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c 




<T 



a 



X 




4-J 




X 



JL 



xy 




a) 

£ 

<D 




P/Dt 



P/Dt 



P/Dt 





X 

R 

t 

FE 


CV 

FE 


CV 

FE 


CV 


TD 

• 

36 

0.100 

0.007 

-0.637 


-0.599 

1.76 


1.76 

-0.792x10 

i-4 

-0.003 


O 

72 

0.220 

0.020 

-0.523 


-0.520 

1.57 


1.58 

-0.018 


0.016 


-a 

o 

108 

0.420 

0.040 

-0.312 


-0.333 

1.11 


1.13 

-0.033 


0.038 


4-> t- 
0) 

144 

0.647 

0.053 

-0.154 


-0.149 

0.670 


0.659 

-0.040 


0.036 


• n 

CLf- 

180 

0.820 

0.073 

-0.056 


-0.089 

0.405 


0.408 

-0.031 


0.031 


L_ 4— 

a> 

CL 

216 

0.947 

0.080 

-0.010 


-0.243 

0.271 


0.259 

-0.024 


-0.018 

ro 

T3 













CO 

03 

O 

28 

0.007 

0.100 

-0.618 


-0.624 

1.84 


1.84 

-0.004 


-0.003 

CTi 

, — • 

64 

0.020 

0.220 

-0.657 


-0.638 

1.92 


1.94 

-0.015 


-0.023 


O *r- 
TS 

100 

0.040 

0.420 

-0.708 


-0.680 

2.28 


2.30 

-0.118 


-0.116 


X 

136 

0.053 

0.647 

-0.656 


-0.731 

3.26 


3.35 

-0.555 


-0.459 


a> *r~ 
e 

172 

0.073 

0.820 

-0.099 


-0.185 

4.88 


4.81 

-1.82 


-1.91 


r— 4-5 
fT? 03 
£ 

708 

0.080 

0.947 

4.45 


5.00 

6.65 


3.16 

-3.62 


-5.73 
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Appendix C 

Determination of Stresses Using a Finite-Difference 
Representation of Equilibrium Equations 

In an x-y coordinate system, the plane stress equilibrium 
equations, in the absence of body forces, are: 

(C-l) 
(C-2) 

The stress-optic law for orthotropic photoelastic materials can be 
written as 

+ cf = N t co$ ( 20 t) ' N R cos(2G R ), (C-3) 

- °v T vv N T N R 

— I L = j- sin(29 T ) - sin(20 R ). (C-4) 

The constants Ci - C 5 are calibration constants for the material, Ny 
and Gy are the total isochromatic fringe number and the optical isocline 
angle, respectively, and Nr and e R are the residual isochromatic fringe 
number and the residual optical isocline angle. The latter two 
parameters can be thought of as calibration constants for the material. 

Consider a rectangular array of points in an x-y coordinate system 
as in fig. C-l. For point i , j + point a, the following 
approximations can be written: 


^ ^ V V 

— - + - *2L = o, 
5x ay 


^y V da v 

— + — X- = o 
ax ay 
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(C-5) 


6X 1.JH1 

dT xy 

ay 1 .j4 


= ax ^°x^i,j+l ‘ 


5t 


xy 


1 (!!sx 

2 1 ^ i,j+i ay i,j 


}. 


(C-6) 


where 


and 






3y i.j« 


5-r 


_xy 


ay 


i*J 


?iy t( T xy >i , j +1 ' (l yy j('l, i+l ^ 

<c -» l 


Equation C-l now can be approximated, at point a, as: 


- K^i.j + W {{x xyWj ‘ ^xyVl.j + (c * 9) 

KyVl.j+l " Ky)i-l,j+l} = °* 

Equation C-9 is an approximation to the equilibrium condition in the x 
direction. It is applied at specific points along the midrow of the 
three-line mesh, row i. This expression averages shear stresses. 

For point (i + poi nt b, the following approximations may be 

written : 


where 


and 


do y 

ay 

^ CT y ^i+i.j * *- CT y ^ 

(C-10) 

!Z5* i . 

dx i+ 2 *»J 

= i + a ^y }, 

2 ,ax i + u ax i,j 

(C-ll) 

® T xy 

ex i+l,j 

2ax ^ T xy^i+l,j+l ^ T xy ^i+l ,j-l ^ 

(C-l?) 

a-c xy 
ax i,j 

= 22sx ^xy^i.j+l * ^ T xy^i ,j-l^’ 

(C— 13) 


Equation C-2 can be approximated at point b as 
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(C— 14) 


( a y)i+l,j ’ (°y^i,j + 4ax ^ T xy^i+l,j+l " ( T xy)i+l,j-l 
+ ^ T xy )i ,j+l " ( T xy )i ,j-l J = °* 

Equation C-14 can be applied at discrete points b and b' along the 
columns of the mesh, columns ...j-2, j-l,...,j+2... The equation also 
averages shear stresses. 

Equations C-3 and C-4 can be written as follows to put them into 
context: 


(°x)i ,j _ + ^ T xy ji , j = 

c x c 2 c 3 


N T 005(20^- ) - N R cos(20 R ) 

1 9 J 1 »J 


(C-15) 


+ ^ T x,y^i ,j 


'T 


-pJ-sin(2e ) 


- s1n(2e 0 ). 


(C-16) 


(Here N R and 0 R are assumed to be the same at each mesh point). 

Equations C-9, C-14, C-15, and C-16 are the four equations which 
can be used to compute the stresses <j v , a u , and x at each point i,j. 

* j Ajr 

It is obvious that x direction equilibrium is not enforced at the same 
points as y direction equilibrium and here y direction equilibrium is 
enforced at twice as many points. There is nothing special about this, 
it is simply a direct result of the mesh geometry, i.e. long and thin. 
The photoelastic equations should be applied at every point. The equil- 
ibrium equations can be applied so there are exactly as many equations 
as there are unknown stresses. Alternatively, the equilibrium equations 
can be applied frequently enough so that there are more equations than 
there are unknown stresses. Solving an overdetermined set of equations 
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will give stresses that satisfy all four equations in a least-squares 
sense. 

Finally, of course, there are boundary conditions. At some points 
one or more stresses may be known, e.q. 5 = 0. These are simply 

more equations in the equation count and can be used to give a unique 
solution or an overdetermi ned solution. 
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